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abstract 


-Three  models  are  examined  to  study  the  transmission  of  ocean  waves 
through  an  ice-field.  In  each  case  the  effect  of  ice  thickness,  water 
depth,  and  the  wavelength  and  angle  of  incidence  of  thi  incoming  ocean  wav 
is  considered.  In  Model  1  the  ice  is  assumed  to  consist  of  floating  non¬ 
interacting  mass  elements  of  varying  thickness  and  the  sha! lew-water 
approximation  is  utilized  to  simplify  the  equations.  A  simple  cosine 
distribution  varying  in  one  direction  only  is  assumed.  In  Model  !l  '■ 
mass  elements,  of  constant  thickness,  interact  through  a  bending  stiffness 
force  so  that  the  ice  acts  as  a  thin  elastic  plate.  The  mass  elements  are 
connected  through  a  surface  tension  force  in  Model  111  so  that  the  Ice  is 
simulated  by  a  stretched  membrane.  In  both  Models  II  arc?  iil  the  full 
linearized  equations  are  solved.  Because  of  the  complexity  of  the  res  sit¬ 
ing  analysis,  calculations  of  the  reflection  and  transmission  coefficients 
and  the  pressure  under  the  ice,  are  made  in  Model  t|  on  the  basis  of  the 
shallow  water  approximation. 
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NOMENCLATURE 

Unless  indicated,  equation  numbers  refer  to  Part  II  of  the  report. 

A  constant  defined  by  Eq.  (3-5);  constant  defined  by  Eq .  (5.22) 

A.(i  =  0,1,2)  constants  defined  by  Eq.  (5.25) 

A.  .(i,j  =  '.2)  constants  defined  by  Eq.  (3-56) 

J 

a  constant  in  Mathieu's  equation,  defined  by  Eqs.  (2.7) 

and  (3-2)  of  Part  I 

—  iaQ  pure  imaginary  roots  of  Eq .  (2.1) 

±a  (n  =  1,2,3...)  real  roots  of  Eq.  (2.1) 
n 

±a  ’  -  Ha*  -  k2)* 

o  o 

B  constant  defined  by  Eq .  (2.2) 

±ibo  pure  imaginary  roots  of  Eq .  (2.2) 

±bn(n  =  1,2,3...)  real  roots  of  Eq.  (2.2) 

±b  1  =  ±(b  3  -  k2)^ 

o  o 

C  path  of  integration  for  the  integral  in  Eg.  (3-33) 

C',  C1 1  deformed  paths  of  integration 

Cx  ,  0%  constants  defined  by  Eq.  (3.30) 

c  =  -^gH  shallow  water  wave  velocity 

c.(i-  1,  2,3,4)  constants  defined  by  Eq.  (3-19) 
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±c 


o  ’ 


±c 

o 


±c  1 

o 

±c  ' 

o 


D 


D 

E 


M<*) 


F(x,y) 

G(y) 

g 

g(P) 

H 

H(x) 

h 


complex  roots  of  Eq.  (2.2) 


-  ±  (c  3  +  ka)2 


±  (S,3  + 


Eh 


12(1  -v3  ) 


;  strip  in  the  complex  or-plane, 


region  t  >  -  k  in  the  complex  a-plane 
region  t  <  0  in  the  complex  <*-plane 

Young's  modulus  for  ice 
iAK  (-  b  *) 


=  J(«)  - 


+ _ o 


Or+b 


o 

defined  by  Eq.  (1.8)  of  Fart  I 
defined  by  Eq .  (2.6)  of  Part  I 
acceleration  due  to  gravity 
defined  after  Eq.  (3-53) 


k  <  T  <  0 


depth  of  water 

solution  of  Eq.  (1.10)  obtained  from  separation 
of  variables 

constant  ice  thickness 

mean  and  fluctuating  ice  thicknesses  in  Part  1 


variable  ice  thickness  in  Part  1 


v  i  i  i 


h(x,y) 
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I  coefficient  of  incident  potential 

i  =  \J-  1 

J{Qf)  =  cia  +  C2 

K  wave  number  for  incident  waves  in  deep  water 


K_  wave  number  for  incident  waves  in  shallow  water 

V 


=  (K  S  -  k2 


K±(a)  functions  regular  and  non-zero  in  respectively, 

”  defined  by  Eq.  (3-25)  and  evaluated  in  Appendix  B 

I  wave  number  for  undulations  in  the  ice  (Part  l) 

|  wave  number  for  x-component  of  incident  waves 

k  (n  =  0,1,... 5)  roots  of  the  equation  Mk  6  +  (i-Ljk.,3  =  K  3 
n  n  11  o 

V  -  (C  -  ^ 

p.u^h 

I  * 


M 

P 

<P 

P 

po 

q 


D/pg 

coefficient  in  pressure  fluctuation  term,  Eq.  (5.38) 

non-d  imer.s  Iona  1  value  of  the  pressure  amplitude 
water  pressure 
atmospheric  pressure 

constant  in  Mathieu's  equation  defined  by  Eq.  (2.8) 
of  Part  I 
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R 


s 

s 

T 

f 

T 

s 

t 


*,y 


z 


coefficient  of  reflected  potential 
reflection  coefficient 
=  ly3  +  (z-H)3]2 
=  T./pg 

=  p./p  ,  specific  gravity  of  ice 
coefficient  of  transmitted  potential 

transmission  coefficient 
suriace  tension  force 

t  ime 

horizontal  space  co-ordinates 
vertical  space  co-ordinate 


a  =  a  +  iT  ,  complex  Fourier  transform  variable 

@  constant  defined  by  £q.  (3.10) 

£«  =  (f  -  k2)^ 

V  =  (a2  +  k2)^ 

e  _  1  +JI' 

2 

fl  defined  by  Eq .  (2.6) 

G  angle  between  direction  of  incident  wave  and  normal  to 

ice  field 

3.J.  angle  between  direction  of  transmitted  wave  and  normal 

to  ice  field 


x 
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e  . 

cr  it 

\ 

k  . 

crit 

k  . 

min 


v 


P 


a 


®(*5y,t),  0(y,z) 


Y  ,  Y 
+  * 


Hy.z) 

U> 


critical  angle  at  which  complete  reflection  occurs 
wavelength  of  incident  wave 

critical  wavelength  above  which  all  waves  penetrate 
the  fee 

minimum  wavelength  below  which  effectively  no  penetration 
occurs  (Part  I) 

Poisson's  ratio  for  ice 

time  dependent  surface  elevation 

dens ity  of  water 

density  of  ice 

real  part  cf  complex  variable  a 
imaginary  part  of  complex  variable  a 
time  dependent  velocity  potential 

shallow  water  velocity  potential  defined  by  Eq .  (1.5) 
in  Part  I;  defined  by  Eq.  (1.9)  in  Part  II,  respectively 

Fourier  transform  of  i>(y,z) 

half-range  Fourier  transforms  defined  by  Eqs.  (3.1^) 
and  (3.15) 

ib  'y 

=  ¥(y.z)  -  e  0  cosh  b  z 


wave  frequency 
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DIVISION  OF  STUDY 

The  study  is  divided  into  two  distinct  parts.  Part  l  consists  of 
Model  I,  the  propagation  of  waves  through  ice  of  variable  thickness  in 
shallow  water.  This  part  is  complete  in  itself,  containing  the  analysis, 
some  numerical  results,  and  a  qualitative  discussion  of  the  importance  of 
the  various  parameters  in  predicting  wave  transmission  and  reflection. 

The  main  part  of  this  work  was  carried  out  by  Professor  T.  V.  Davies, 
Visiting  Scientist,  Davidson  Laboratory,  September  1965  to  July  1966. 

Part  il  contains  Models  II  and  lli,  where  the  ice  is  assumed  to  have 
a  flexural  stiffness  and  a  surface  tension-type  force,  respectively.  For 
Model  II,  numerical  results  for  transmission  and  reflection  coefficients, 
together  with  the  amplitude  of  the  pressure  fluctuation  on  the  bottom 
under  the  ice,  are  given,  together  with  a  discussion  of  the  results. 

Each  part  is  complete  in  itself  and  may  be  read  separately. 
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INTRODUCTION 


This  report  forms  the  theoretical  part  of  a  combined  theoretical 
and  experimental  study^  into  the  effect  of  water  waves  on  an  ice-field  in 
water  of  finite  depth. 


Wave  transmission  and  reflection  in  finite  and  infinite  depths  of 
water,  partially  ice-covered,  have  been  the  subject  of  a  number  of  theoret¬ 
ical  stud ies' ,6  contrast  to  the  scarcity  of  experimental  ones.^’^ 
The  theoretical  studies  have  been  dominated  by  the  basic  assumption  that 
the  sheet  of  ice  can  be  represented  by  a  semi-infinite  rigid  sheet,  or  by 
a  sheet  composed  of  non-interacting  floating  point  masses.  So  far,  theory 
reveals  that  the  transmission  of  propagating  undamped  waves  under  the  ice 
depends  upon  the  assumed  surface  condition,  the  ice  thickness,  the  angle 
of  incidence  of  the  incident  waves,  the  wavelength  of  such  waves,  and  the 
depth  of  water. 


Heins  has  assumed  the  ice  to  be  a  continuous  rigid  sheet  without 

movement,  extended  over  a  semi- inf inite  region  of  finite  depth;  he  studied 

2 

the  water  wave-ice  interaction  under  these  conditions.  Peters  assumed 
that  the  'ce  consisted  of  broken  pieces  with  no  interaction,  i.e.,  neither 
stiffness  nor  elasticity,  extending  over  a  semi-infinite  region  on  the 
surface  of  an  infinite  depth  of  water;  he  found  that  there  is  a  critical 
incidence  frequency  wh<:t  determines  whether  or  not  there  will  be  a  damped 
transmitted  wave  srier  the  interaction.  Keller  andWeitz  have  analyzed 
the  same  problem,  W.  with  water  of  finite  depth.  Shapiro  and  Simpson 
have  made  a  numerical  study  of  the  above  reference  which  shows  that  water 
waves  entering  an  ice  field  are  damped  exponentially  with  increasing  ice 
thickness.  Keller  and  Goldstein"*  have  considered  the  reflect  ..,i  of  water 
waves  from  a  region  covered  by  floating  matter  with  and  without  surface 
tension;  they  have  solved  the  wave-ics  interaction  problem  for  an  arbitrary 
incidence  angle  by  utilizing  the  shallow  water  theory.  Their  study  reveals 
the  importance  of  the  incidence  angle  as  well  as  of  the  degree  of  stiffness 
of  the  floating  materia).  Perhaps  the  most  realistic  model  for  the 
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transmission  of  water  waves  by  large  floes  Is  given  by  Stoker  although 
his  motivation  was  different.  He  was  concerned  with  the  effectiveness  of 
a  floating  beam  having  a  known  flexural  stiffness,  as  a  breakwater.  The 
treatment  is  restricted  to  two  dimensions  in  the  sense  that  the  incident 
waves  enter  the  ice  normally,  and  the  shallow  water  approximation  is 
ut i 1 ized. 

7  8 

Finally,  the  observations  of  Robin'*  on  wave  propagation  through 
ice  fields  are  considered  to  be  fundamental  in  which  he  demonstrates  the 
existence  of  a  higher  cut-off  wavelength  depending  on  the  size  of  floes, 
above  which  no  attenuation  of  water  waves  occurs,  and  a  lower  cut-off 
wavelength  below  which  the  incident  wave  is  completely  absorbed. 

The  present  study  consists  of  three  distinct  models.  In  the  first 
model  the  ice  is  assumed  to  be  made  up  of  floating  non-interacting  elements 
of  varying  mass  density  distribution.  The  shallow  water  approximation  is 
used  and  a  simple  thickness  distribution  is  considered  to  simplify  the 
analysis.  A  qualitative  discussion  is  made  of  the  Importance  of  the  varia¬ 
tion  in  floe  thickness,  the  angle  at  which  the  incident  wave  approaches 
the  ice-field,  the  wavelength  of  the  incident  wave,  and  the  depth  of  water. 
It  is  shown  that  transmission  of  a  g'ven  incident  wave  at  a  given  incidence 
angle  is  critically  dependent  upon  the  parameters  of  the  problem,  but  that 
waves  which  are  long  enough  will  always  be  transmitted  into  the  ice.  In 
addition,  some  computations  are  made  which  indicate  in  effect  a  wavelength 
below  which  no  incident  waves  penetrate  the  ice.  These  results  confirm  the 
observations  made  by  Rob  in.  ^ 

In  Model  It  the  i-ce  rs  assumed  to  have  a  constant  thickness  h  ,  but 
is  allowed  to  bend  tc  permit  the  transmission  of  water  waves.  This  corres¬ 
ponds  to  Stoker's  niodel,  but  in  this  case  waves  may  be  incident  on  the  ice 
from  any  angle  and  the  problem  is  solved  for  any  depth  of  water  H  .  Be¬ 
cause  of  the  complexity  of  the  analysis,  numerical  computations  are  made 
only  under  the  assumption  of  the  shallow  water  theory.  It  is  shown  that 
a  critical  incidence  angle  exists  for  each  wavelength  above  which  the 
incident  wave  is  completely  reflected  by  the  ice-field,  without  transmission 
taking  place.  In  addition  there  exists  a  critical  wavelength  above  which 
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all  incident  waves,  at  any  incident  angle,  penetrate  the  ice-field  as  un¬ 
attenuated  transmitted  waves  of  reduced  amplitude.  Extensive  calculations 
are  made  which  emphasize  the  importance  of  wavelength,  incidence  angle, 
ice  thickness  and  water  depth,  in  determining  transmission  and  reflection 
coefficients,  and  the  pressure  fluctuations  on  the  bottom  under  the  ice. 

It  is  felt  that  this  model  provides  a  good  representation  of  the  trans¬ 
mission  of  waves  through  large  ice  floes  which,  according  to  Robin, ^  have 
been  observed  to  bend  to  allow  waves  to  propagate  through. 

In  Model  li!  the  ice  is  assumed  to  be  made-up  of  floating  point 
masses  which  are  connected  through  a  surface  tension  force.  This  model 
was,  in  fact,  considered  before  the  more  complicated  Model  II,  and  many 
of  the  difficulties  which  arise  in  this  more  realistic  case  were  first 
encountered  and  solved  in  the  simpler  boundary-value  problem  of  Model  III. 

It  is  acknowledged  that  this  surface  tension  model  is  not  a  realistic 

representation  of  the  wave-ice  interaction  problem  but  it  is  of  some 

5 

academic  interest  since  it  extends  the  work  of  Keller  and  Goldstein  to 
finite  depth  of  water,  it  is  only  given  h-iaf  treatment  since  it  follows 
closely  the  techniques  used  in  Model  II. 
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PART  I 

MODEL  l.  WAVE  PROPAGATION  THROUGH  AN  ICE-FIELD 
OF  VARYING  THICKNESS 

!.  INTRODUCTION  AND  EQUATIONS  GOVERNING  THE  PROBLEM 

A  semi-infinite  sheet  of  Ice  of  variable  thickness  is  at  the  surface 
of  an  ocean  of  uniform  depth  H  ;  the  sheet  occupies  the  domain  o  <  y  <  , 

-  n'  <  x  <■  ,x>  ,  z  =  H  and  the  remainder  of  the  domain  z  =  H  is  open  water. 
Ocean  waves  are  assumed  to  impinge  normally  on  the  edge  of  the  ice-sheet 
and  the  problem  arises  of  the  nature  of  the  transmitted  wave.  Here  we 
assume  from  Keller  that  the  problem  can  be  approached  using  the  approxi¬ 
mations  of  shallow  water  theory,  that  is,  we  are  assuming  that  the  wavelength 
^  of  the  incident  ocean  waves  is  long  compared  with  depth  H  of  water. 

In  the  first  place,  if  we  approach  the  problem  on  the  basis  of  the  lin¬ 
earized  theory  only,  the  velocity  potential  $(x,y,z,t)  for  the  liquid 
motion  satisfies 

4  +4  +4  =  0  ,  0  <  z  <  H  (1.1) 

xx  yy  zz 

and 

4z  =  C  ,  x  -  C  (1.2) 

At  the  free  surface  of  the  ocean  where  §(x,y,t)  is  the  elevation  above 
the  undisturbed  level,  we  have,  on  the  basis  of  linearized  theory, 

d£  _  34  _  _  34 

dt  "  “  dz  ’  9*  “  dt 

so  that 

4^  +9^  =  0  ,  z=H,-«<y<0  (1.3) 
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The  difference  of  pressure  across  the  ice-sheet  is  given  by 

P  -  -  P(||  ‘  9^ 

and  the  equation  of  motion  of  an  ice  element,  in  the  absence  of  bending 
stiffness  and  surface  stress  effects,  is 

p  -  Po  =  pth(x,y)Stt 

where  p.  is  the  ice-density  and  h(x,y)  the  variable  thickness  of  the 
ice.  Hence,  the  condition  to  be  satisfied  on  the  ice-sheet  will  be 

H  -  =  s  h(x>yKtt 

(1-4) 

a?  =  _  aa 

5"t  dz 

where  s  =  p./p  is  the  specific  gravity  of  ice. 

In  order  to  simplify  the  problem,  we  now  invoke  the  shallow  water 
theory  approximations  and  we  do  so  by  writing  (see  Keller"*) 

£(x,y,z,t)  =  <p(x,y,t)  -  iz3^  +  <Pyy)  +  0(z)3  (1.5) 

This  expression  will  satisfy  Eq.  (l.l)  to  the  first  order  and  will  also 
satisfy  Eq .  (1.2).  In  order  to  satisfy  Eq.  (1.3),  the  function  v(x,y,t) 
must  satisfy 

9tt  -  9  H  (cp^  +  <pyy)  =  0  ,  -®  <  y  <  0  (1.6) 

and,  after  some  reduction  of  Eq.  (1.4),  we  find  that  9(x,y,t)  must 
sat isfy 

<?tt-  9  ^(4^+  9yy)-  s  Hh  (x,y)(|pr  +  §^3-}<Ptt=  0,  0  <  y  <  ®  0-7) 
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Equations  (1.6)  and  (1.7)  constitute  the  basic,  equations  of  the 
problem,  the  function  h(x,y)  being  the  prescribed  variable  thicknes:  of 
the  ice. 

We  can  take  the  function  t?  to  contain  the  time  through  n  exponen¬ 
tial  factor,  and  if  we  write 

cp®  F(x,y)  exp(-iwt)  (1.8) 


then  we  have 


F  +  F  +  —  F  =  0  _eo  <  y  <  Q 
xx  yy  gH  *  1 


0.9) 


and 


'  2  >  2 

•j  1  -  — —  h(x,y)r  (F  +  F  )+^F  =  0  ,0<y<®  (1.10) 

k  g  ’"J  '  xx  yy  gH 


At  the  boundary  y  =  0  ,  it  will  be  necessary  to  have  9  and  9^  contin¬ 
uous,  expressing  the  continuity  of  normal  and  transverse  velocity,  hence 


F  . 


F 

y 


continuous  at 


y  =  0 


(l.li  i 


2.  ICE  THICKNESS  VARYING  IN  y-DIRECT[0N 
WITH  NORMALLY  INCIDENT  WAVES 

Here  we  take  the  ice  thickness  h  to  vary  only  in  the  y-direction 
and,  in  order  to  simulate  the  effect  of  ice-floes,  it  will  be  assumed  to 
have  the  following  structure: 

h  =  +  hx  cosky  (2.1) 

where  h,  s  h  and  2ti/k  is  the  wavelenqth  of  undulations  in  the  ice. 
i  o 

Ue  look  for  solutions  in  which  the  incident  ocean  waves  impinge 
normally  against  the  edge  of  the  ice  sheet,  in  which  case  the  complete 
expression  for  9  in  -  co  <  y  <  0  can  be  taken  to  be 
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cp  =  exp 


\iuJ(c  -  *>} 


(2.2) 


where 


The  transmitted  wave  in  C  <  y  <  “  will  be  of  the  form 


cp  =  G(y)  expi-)<*>t) 


(2.3) 


and  from  (1.10)  it  follows  that  G(y)  satisfies  the  differential  equation 


sh  UJ3 
o 

9 


■)  - 


sl^ui3 

g 


cos 


0 


(2.4) 


It  is  not  necessary  to  work  with  this  differential  equation  in  the  above 
exact  form  since,  when  we  insert  typical  value  of  the  constants  h^  ,  ht  , 
s  ,  ou  ,  we  find  that 


sh,  u.2  sh  (B3 

— —  «  l  ,  — —  «  1 

9  S 

for  waves  whose  periodic  time  is  greater  than  4  seconds  and  we  shall 
assume  that  the  investigation  is  restricted  to  this  class  of  waves.  In 
this  case  it  is  permissible  to  expand 


{< 


sh  uj3 


1  - 


-) 


shjUf3 


kyj- 


in  the  form  of  a  convergent  infinite  series  in  the  parameter  sh ^/(g-sh^w2 ) 
and  thus  to  write  (2.4)  in  the  approximate  form 


u? 

_2H 


fi¬ 


sh 

-4-> 


{> 


sh  jU) 


0- 


cos 


0 


(2.5) 
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Accordingly,  if  we  now  write 


ky  =  2T| 


(2.6) 


Equation  (2.5)  becomes 

j2  r 

+  (a  +  2q  cos  2H)G  =  0 


(2.7) 


(2.8) 


(2.9) 


Q 

which  is  the  standard  form  of  the  Mathieu  differential  equation  (McLachlan, 
Theory  and  Application  of  Mathieu  Functions).  Periodic  solutions  of 
Mathieu's  equation  e.-ist  only  under  special  circumstances  and  in  order  to 
make  this  clear  we  refer  to  the  stability  diagram  in  which  q  is  the 
abscissa  and  a  the  ordinate  (Fig.  I). 

If  the  values  of  q  and  a  are  such  that  the  representative  point 
(q,a)  in  the  stability  diagram  lies  on  the  curves  denoted  by 


s 

o 


b 

i 


there  will  be  a  periodic  solution  of  Eq.  (2.9).  When  q 
small,  it  is  known  (see  McLachlan)  that  the  equations  of 
are  as  follows: 


i  s 
a 


o 


sufficient 

'  ’  ai 


a 

o 


+  I2S  q4+  °(q6) 


(2.10) 
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q "  i q2+  w  q3_  rhs  q4+  0(qS) 

(2.11) 

q  "  S’  q3"  s?  q3"  755^  q4+  0(q6) 

(2.12) 

12  a  +735?Tq4  +  0(qS) 

(2.13) 

Tr5’  -w;<*  +  °W“> 

(2.14) 

and  the  solutions  corresponding  to  ac  ,  ^  ,  ...  are  as  follows 

in  the  usual  notation 


s  ; 

0 

CO 

ceo(T|)  -  2 

0 

A(0) 

A2n 

cos  2nfl 

(2  15) 

•v 

se.(Tl)  =  2 
o 

r0  ) 

2n+l 

sin(2n+l)T] 

(2.16) 

a  ; 

i 

ce  (11)  =  S 

!  o 

A0) 

Zn+I 

cos  (2n+l)T) 

(2.17) 

Elsewhere  in  the  stability  diagram  the  solution  for  G(Tj)  of  Eq.  (2.9) 
is  of  the  form 

G(T,)  =  e1^  $  (Tl.cr)  +  e-TltA(a)  4  (T] , -cr)  (2.18) 

where  $(7,11)  is  a  periodic  function  of  11  ,  is  the  characteristic 

exponent  of  the  Floquet-Poi.ncare  theory,  and  cr  is  a  parameter.  We  can 
illustrate  the  significance  of  the  parameter  <T  by  considering  first  of 
all,  the  area  of  the  stability  diagram  lying  between  a^  and  bg  . 
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Whittaker  obtained  the  form  of  solution  in  this  domain  and  he  finds  that 
(McLachlan,  p.  70)  one  solution  of  (2.9)  is  given  by  e”^  ^  as 

fol lows: 


$(z,cr)=  sin(z-o)+  s,sin(Bz-o)+  c  cos(Jz-a)+  sin (5z.*<J)+* * •  (2.19) 

o  3  5 

and  a  and  n  are  expressed  in  terms  of  q  and  cr  as  follows: 

a(o,q)  =  l-q  cos2cr+ q2  (-1+ icos4o)+ gjj- q3cos2o  +  ...  (2.20) 


M.{cr,q)  =  -  y  q  sin2cr  +  -~g  q3sin2a  +  ..  .  (2.21) 

li  will  be  noted  that  a(0,q)  -  bj  ,  a(-  -z  ,  q)  =  at  ,  where  bx  ,  a,  are 
the  values  given  in  (2.11)  and  (2.12).  Whittaker  finds  that  the  value 

cj  =  -  y  n  +  j0  ,  6  real  ,  ^  0  (2.22) 

gives  the  region  between  ax  and  and  here  H  is  purely  imaginary. 

This  is  therefore  designated  a  stable  area  in  the  stability  diagram. 

In  a  similar  way,  the  value 

a  =  ie  ,  S  real  ,  *  0  (2.23) 

qives  the  stable  region  between  a  and  b,  .  The  area  between  b  and 

O  l  1 

is  such  that  a  is  real  and  lies  in  the  range  , 

-  ~  s  a  si  0  (2.24) 

and  this  is  designated  an  unstable  area  in  the  stability  diagram.  The 
remainder  of  the  Stability  diagram  can  be  investigated  and  described  in  a 
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s iml lar  way . 

The  solution  given  in  Eq .  (2.18)  indicates  that  the  second  solution 
can  be  deduced  from  the  above  described  solution  merely  by  a  change  in  the 
sign  of  o  . 

Returning  now  to  the  wave  problem,  we  see  that  as  long  as  the 
representative  point  (q,a)  lies  in  any  unstable  domain  (shown  shaded  in 
the  stability  diagram),  the  corresponding  solution  for  G(f|)  must  be  an 
attenuated  transmitted  wave;  the  degree  of  attenuation  will  depend  upon 
the  magnitude  of  p.(°0  and  we  see  that  in  Eq.  (2.21),  p  will  attain  its 
maximum  value  in  the  unstable  domain  (bt  ,  ^  )  for  cr  =  -  ^  provided  q 
is  sufficiently  small.  On  the  other  hand,  if  the  representative  point 

(q,a)  ties  on  ,  _  or  within  the  stable  unshaded 

areas  (a^  ,  b^)  ,  ( a L  ,  bg)  ,  ...  ,  the  solution  for  G( T))  v.'ill  be  an 
undamped  transmitted  wave.  The  character  of  this  undamped  transmitted: 
wave  will  be  periodic  (but  not  a  sine  wave)  when  (q,a)  actually  lies  on 
the  curves  a  ,b  ,a  ,b  the  structure  when  (q,a)  lies  in  a 

Stable  area  Such  as  (aQ  ,  b^)  is  not  periodic  but  almost  periodic.  This 
description  of  regions  of  damped  and  undamped  solutions  can  be  expressed 
in  a  most  useful  form  by  returning  to  the  definitions  of  a  and  q  in 
Eqs.  (2.7)  and  (2.8).  It  will  be  noted  that  each  of  these  parameters 
contains  the  geometrical  quantities  H  ,  h  ,  h ^  ,  the  specific  gravity 
s  as  well  as  the  wavelength  7^/k  of  the  ice  undulations  and  the  ocean 
wave  freouency  parameter  a,-;  itwill  be  observed  also  that  the  parameter  cu  , 
when  eliminated  between  ;2.7)  and  (2.8),  leads  to  the  relation 


=  (■ 


sk  hj^H' 


(2.25) 


For  a  fixed  value  of  (S/s^hjH)  ,  this  relation  will  be  represented  by  a 
parabola  I’  .  This  parabola  will  intersect  the  curves  b^  ,  a^  ,  b^  ,  a^  , 
in  points  8  ,  ,  B g  ,  Ag  ,  . . .  as  shown  in  the  diagram,  and  it  is  a 

straight-forward  matter  to  locate  these  points  analytically  when  q  is 
not  large.  It  is  then  possible  to  determine  the  values  of  at  the 
positions  ,  A^  ,  B  ,  A^  ,  ...  which  we  can  denote  by  ^(B^)  ,  ^(A^)  , 
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It  then  follows  that  the  domains  of  w  defined  by  0  <  tu  <  )  , 

tu (A x )  <  a)  <  cu(Be)  ,  ...  lead  to  stable  or  undamped  transmitted  waves, 
while  the  ranges 


iu^)  <  u»  <  to(Al)  ,  is(Bs)  <  «>  <  ao(Aa)  ,  ... 


give  damped  or  attenuated  transmitted  waves.  The  breadth  of  these  ranges 
will  depend  upon  the  parameter  (8/sksh  H)  which  is  seen  therefore  to  play 
a  crucial  role.  If  we  write 

r  =  — ■— —  (2.26) 

sk  h  H 

i 

the  parabola  T  has  the  equation  a*  =  2rq  and  the  intersection  of  T 
with  the  curve  b  ,  for  example,  can  be  obtained  by  combining  (2.26)  and 
(2.11),  namely 


a  =  I  -  q  +  0(qs ) 


(2.27) 


so  that  the  intersection  is  given  by 

2rq  =  (I  -  q)2 
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that  is 

q  =  (1  +  r)  -  yj  r(2  +  r) 

and  the  value  of  tu  then  follows  using  (2.8).  The  approximate  positions 
of  cu( B1 )  ,  ^(A^,  ...  can  also  be  determined  in  this  way. 

To  summarize  the  conclusions  of  this  section  of  Part  I,  the  impor¬ 
tant  feature  of  the  results  is  that  when  all  the  geometrical  parameters 
are  fixed,  namely,  the  depth  of  water,  average  ice  thickness  and  departure 
of  ice  thickness  from  the  mean,  the  wavelength  of  undulations  in  the  ice, 
there  can  be  certain  wavelength  ranges  for  the  incident  ocean  waves  which 
will  be  attenuated  on  the  ice-sheet  and  other  wavelength  ranges  which  will 
be  transmitted  as  undamped  waves. 


3.  MODIFICATION  FOi  OBLIQUELY  INCIDENT  VAVES 

The  above  analysis  requires  modification  whenever  the  incident  wave 
impinges  upon  the  ice-field  at  an  angle  other  than  zero.  In  order  to 
carry  out  this  modif ication,  we  refer  to  Eqs.  (1.9)  and  (1.10).  Clearly 
the  solution  of  Eq.  (1.9)  representing  an  incident  wave  making  an  angle  9 
with  the  ice-field  is 

i K  y  cos  8  +  iK  x  sin  6 
„  o'  o 


where 

k  -  _HL_  -  2H 

Ko  ~  VgH  ~  x 

and  X  is  the  wavelength  of  the  incident  wave. 

Assume  a  solution  of  Eq .  ( 1 _ 1 0)  of  the  form  H(x)  G(y)  . 
stitution  into  (1.1 0)  gives 


h"i 

1  _  (g”< 

ft 

l — 

b 

7* 

HI 

'  1 

l G< 

ft 

1  1-SHKoah(y)j 
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where  the 
Since  the 
side  is  a 
To  obtain 
necessary 
solutions 
x  ,  -  “  < 
y  >  0  as 


thickness  of  the  ice  is  assumed  to  vary  in  the  y-direction  only, 
left-hand  side  is  a  function  of  x  only,  while  the  right-hand 
function  of  y  only,  each  side  must  be  equal  to  a  constant, 
an  oscillatory  x-variation  for  all  x  ,  -  <  x  <  ,  it  is 

that  this  constant  be  negative.  Further,  in  order  that  the 
for  y  <  0  and  y  >  0  might  be  continuous  at  y  =  0  for  ai! 
x  <  “  ,  it  is  necessary  for  the  x-dependence  to  be  the  same  for 
for  y  <  0  .  Thus 


H“(x)  = 
H(x) 


.  s 


s  in 


9 


and 


iK  x  sin  8 
H(x)  =  e  ° 


Thus  the  equation  satisfied  by  G(y)  is 


+  K  - - - 

dy  °  Ll-  sHK  2h(y) 
o 


-  sin2e]  G(y) 


0 


(3-D 


Substituting  the  assumed  thickness  distribution  h(y) 
into  Eq.  (3-0  and  making  the  same  approximations  as  before 


=  h  +  h  cos  ky 
o  1 

namely 


sh,  K  aH  <  <  1  ,  sh  K  aH  <  <  1 
1  o  o  o 


gives 


- 

-  +  (a  +  2q  cos  2T|)  G(T|)  =  o 

dip 


where  in  this  case 


sin8  9 


(3.2) 
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2K  3  sh.K  2  H 
o  _ ]_o _ 

k3  (I  -  sh  K  3H)S 
o  o 


(3.3) 


and  the  substitution  icy  =  2T|  has  been  made.  Once  again  Mathieu's  differ¬ 
entia)  equation  is  obtained.  Here  also  the  frequency  dependence  may  be 
eliminated  to  obtain  a  relation  between  a  and  q  ,  but  since,  in  this 
case,  the  relation  is  no  longer  simple,  there  seems  to  be  no  advantage  in 
doing  this. 


4.  DISCUSSION  OF  RESULTS 

(a)  Existence  of  Energy  Cut-Off  And  its 
Dependence  on  the  Various  Parameters 

The  present  mathematical  model,  based  on  the  shallow-water  approxi¬ 
mation,  with  ice  undulations  following  a  cosine  distribution,  clearly  does 
not  represent  an  actual  ice-field.  However,  a  number  of  qualitative  re¬ 
sults  may  be  deduced  which  might  be  expected  to  hold  true  for  more  realistic 
distributions  of  ice  floes.  Thus  it  has  been  shown  that  for  fixed  values 
of  k  ,  h  ,  H  ,  there  may  or  may  not  be  undamped  transmission  of  a  given 
incident  wave  into  the  ice-field.  By  undamped  transmission  is  meant  a 
wave  travelling  through  the  Ice  undiminished  in  amplitude  with  distance  into 
the  ice.  An  attenuated  transmitted  wave  is  a  wave  which  decays  exponen¬ 
tially  and  travels  only  a  short  distance  into  the  ice.  Thus  for  6=0, 
whenever  the  parabola  given  by  Eq .  (2.25)  crosses  a  shaded  region  in 
Fig.  1,  the  transmitted  wave  is  attenuated,  whereas  when  the  parabola 
crosses  an  unshaded  region,  the  wave  will  be  unattenuated  and  will  proceed 
through  the  ice.  Thus  the  domains  of  uj  defined  by 

0  <  uj  <  u^)  ,  t«(Ai)  <  w  <  uu(B3 ) ,  ... 
denote  stable  or  undamped  transmitted  waves,  while  the  ranges 
cu(Bx)  <  (U  <  m(Ax)  ,  tu(B2)  <  w  <  iu(A  J  ,  ... 
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denote  damped  or  attenuated  waves. 

Now  for  shallow  water,  the  relation  between  wavelength  A  and  fre¬ 
quency  UJ  is  given  by 

x  =  2TT<JqH 

(U 

so  that  corresponding  to  the  domains  of  u>  described  above  we  have  the 
wave-hands 

A  >  Mb,)  ,  Max)  >  A  >  \(Ba),  ... 


denoting  wavelengths  of  incident  waves  which  are  transmitted  through  the 
ice  as  undamped  waves,  and  the  wave-bands 

Mb.)  >  a  >  H\)  ,  a(b?)  >  a  >  X(a3),  ... 

corresponding  to  attenuated  transmitted  waves. 

One  striking  observation  which  may  be  made  is  that  there  exists  a 
wavelength 

MB,)  = 

such  that  all  incident  waves  whose  wavelength  satisfies  A  >  A(B  )  are 
transmitted  through  the  ice.  When  A  =  A(B1)  the  first  energy  cut-off 
occurs  and  complete  reflection  of  the  incident  wave  takes  place  for  ail  A  , 
satisfying  A  (A,)  <A<A(B,).  For  A(B2)  <A  <  A(A^)  wave  transmission 
again  occurs.  The  distribution  and  width  of  the  stable  and  unstable 
regions  are  governed  by  the  various  parameters  in  the  problem,  and  it  is 
possible  to  make  the  following  genera!  observat ions . 

As  h  -*  0  ,  corresponding  to  a  uniform  constant  thickness  h  , 

1  o 

then  q  —  0  and  Mathieu‘s  equation  degenerates  into  the  equation 

aG(H)  =  0 


2n"VqH 

~mb7t 
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whose  solution  does  not  of  course  exhibit  stable  and  unstable  regions. 

This  may  also  be  seen  by  considering  the  parabola 

a3=-^Hq  (8=0)  (3-U) 

For  small  hx  the  parabola  becomes  steeper  and  the  width  of  the  unstable 
regions  becomes  smaller,  there  being  no  unstable  regions  in  the  limit 
hj  =  0  (see  Sketch  1,  p.!3  ). 

A  similar  argument  can  be  made  for  the  case  of  small  k  ,  corres¬ 
ponding  to  long  wave  undulations  in  the  ice,  and  small  H  ,  corresponding 
to  very  shallow  water.  In  each  case  the  parabola  (Eq.  [3.4])  is  very 
steep  and  the  width  of  the  unstable  regions  are  small,  vanishing  altogether 
in  the  limit  of  H  ,  k  =  0  .  This  also  follows  from  a  consideration  of 
the  limit  of  the  Eq.  (2.5)  for  small  H  and  k  . 

It  appears  from  Eq.  (3-2)  that  the  effect  of  finite  incidence  angle 
S  is  to  reduce  the  value  of  a  ,  while  q  remains  constant,  thus  moving 
points  on  the  stability  diagram  into  regions  of  greater  instability.  For 
angles  close  to  90°,  where  a  is  close  to  zero,  but  still  positive. 

Sketch  1  indicates  that  the  first  stable  region  0  <  u>  <  ^(E^)  js  wide, 
but  that  subsequent  regions  are  predominantly  unstable.  Thus  the  wave¬ 
length  MBj)  at  which  the  first  energy  cut-off  occurs  is  small,  and 
only  narrow  bands  of  wavelengths  smaller  than  MBj)  penetrate  the  ice. 

Further  information  concerning  the  effect  of  the  various  physical 
parameters  upon  wave  reflection  and  transmission  was  obtained  by  making 
some  computations  which  are  described  in  the  following  section. 

(b)  Numerical  Results 

For  the  case  of  normally  incident  waves  (6  =  0°),  some  computations 
were  made  to  determine  the  breadth  of  the  unstable  regions  correspond ing 
to  complete  reflection  of  the  incident  wave.  Dimensions  corresponding  to 
model  sizes  were  chosen  which  were  then  scaled  using  scale  ratios  of  80:1 
and  200:1.  Thus,  ice  thicknesses  of  hQ  =  =  0.25,  1.00  and  2.00  inches 
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were  considered,  with  wavelengths  of  undulations  in  the  ice  (i.e.,  floe 
length)  of  0.50,  1.00,  and  3.00  feet.  The  water  depth  was  taken  to  be 
0-50  feet  and  the  specific  gravity  of  ice  s  =  0.92  .  Then  from  Eq.  (2.25) 
the  quantity  a  was  determined  for  particular  values  of  q  .  The  equa¬ 
tions  of  the  curves  a  ,  fc>  ,  a  ,  b  up  to  b  were  computed  using  the 
O  1  1  Q  2  ° 

tables  in  Appendix  II  of  McLachlarr  for  q  >  1  ,  together  with  the  series 

3 

expansions  given  by  Eqs.  (2 . I 0)  to  (2.14)  for  q  <  1  (see  McLachlan,  p.  lb- 
17)  •  These  curves  were  then  drawn  and  the  intersection  of  the  parabola 


2 


a 


8 

I  3l  M  ^ 

sk  r^H 


with  each  curve  a  ,  b  ,  ...  b  ,  tabulated.  The  intersection  points 
Ol  € 

give  values  of  q  from  which  the  frequency  and  hence  the  wavelength  of 
the  incident  wave  can  he  determined  using  Eq.  (2.8).  These  points  which 
span  shaded  regions  (Fig.  1)  denoting  unstable  solutions  of  Mathieu's 
equation,  determine  bands  of  wavelengths  corresponding  to  incident  waves 
which  are  completely  reflected  by  the  Ice-field.  The  results  are  given 
in  Table  I.  All  wavelengths  which  lie  outside  the  ranges  indicated  in 
the  table  correspond  to  incident  waves  which  are  transmitted  through  the 
ice  as  undamped  waves.  The  largest  wavelength  given  in  each  case  is  A. ) 
and  ail  waves  having  a  wavelength  A.  >  A^)  propagate  into  the  ice  as 
undamped  waves. 

From  Fig.  I,  as  q  and  hence,  w  s  increases  indefinitely,  the  width 
of  the  stable  regions  intersecting  the  parabola  (2.25)  diminishes  until  a 
value  of  q  is  reached  above  which  the  bands  of  stable  frequencies  are 
indistinguishable  from  points.  Since  the  stability  curves  all  cross  the 
real  axis  for  large  enough  q  (see  McLachlan,^  p.  39),  the  number  of 
intersection  "points"  is  infinite.  Associated  with  this  value  of  q  is  a 
wavelenqth  which  we  denote  by  A  .  .  Thus,  we  define  A  .  ,  as  the  wave- 
length  below  which  the  bands  of  stable  wavelengths  intersected  by  the 
parabola  (Eq.  [2.25]),  are  sensibly  points,  corresponding  to  very  little 
penetration  of  the  ice.  Clearly  this  is  somewhat  arbitrary,  depending  as 
it  does  on  the  accuracy  of  computation,  but  A  .  does  provide  us  with  a 
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useful  bound  below  which  only  discrete  wavelengths  penetrate  the  ice-field. 

Similarly,  for  small  q  ,  and  hence  ,  the  width  of  the  unstable 
regions  intersected  by  the  parabola  (Eq.  [2.25])  diminishes  and  reduces  to 
‘'points."  Also  for  large  floes  (k  small),  or  thin  ice  (hQ,  hx  small)  the 
parabola  is  steeper  and  the  bands  of  instability  given  by  the  intersection 
are  diminished. 

In  the  table,  both  stable  and  unstable  bands  which  are  sensibly 

points  are  omitted.  As  an  illustration  of  the  use  of  the  table,  consider 

the  case  of  water  of  depth  100  feet,  and  ice  of  thickness  4.17  feet.  Then 

for  floe  lengths  of  100  feet,  A  .  =  162  there  being  no  wave  bands  (but 

min 

an  infinity  of  wave  "points")  which  penetrate  the  ice  for  A  <  162  ,  whereas 

if  the  ice  thickness  is  16. 67  feet,  A.  .  has  more  than  doubled  to  344. 

m  1  n 

Also,  note  that  all  wavelengths  penetrate  a  floe  600  feet  long  and  4.17 
feet  thick  (apart  from  an  infinity  of  discrete  wavelengths)  but  if  the 
floe  is  16.67  feet  thick,  there  exists  bands  of  wavelengths  which  fail  to 
penetrate  the  ice. 

(c)  Comparison  of  Results  with  Observation 

7  8 

A  unique  study  of  waves  in  pack  ice  was  made  by  Robin^’  during  a 
voyage  into  the  Weddel I  Sea  aboard  RRS  JOHN  BISCOE  in  1959-60-  He  finds 
that  for  floes  of  around  1.5m  thick  and  40m  or  less  in  diameter, 

"...the  main  energy  cut-off  took  place  when  floe  diameters 
were  about  one-third  of  the  wavelength;  little  loss  of  energy 
occurred  when  floes  were  less  than  one-sixth  of  the  wave¬ 
length  across,  while  no  detectable  penetration  took  place 
when  the  floes  were  half  a  wavelength  or  more  in  diameter." 

From  the  table,  the  closest  comparison  can  be  made  for  floes  of 
4.17-ft  thick  and  100-ft  long.  Then  little  penetration  takes  place  for 

A.  <  162  corresponding  to  a  floe  length  of  five-eights  the  wavelength  or 

7  8 

more  as  compared  to  the  half  wavelength  observed  by  Robinl1  It  is  not 
possible  to  estimate  the  energy  loss  for  long  waves  as  this  requires 
knowledge  of  the  transmitted  wave  amplitude  and  hence  the  full  solution 
of  Mathieu's  equation.  However  the  table  indicates  that  the  first  unstable 
region  corresponding  to  an  energy  cut-off  or  complete  reflection  of  the 
incident  wave  occurs  when  A  drops  to  248  feet  or  when  the  floe  length  is 
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about  two-fifths  of  the  wavelenrtr..  This  compares  favorably  with  the 

7  8 

observed  value  of  one  third  of  a  wavelength  given  by  Robin.  > 

5.  CONCLUSIONS 

It  would  appear  that  the  simplified  model  considered  here  confirms 
qual itatively  the  following  observations  of  Robin: 

(1)  The  existence  of  a  critical  wavelength  at  which  a  major 
energy  cut-off  occurs. 

(2)  The  existence  of  a  wavelength  below  which  very  little 
penetration  of  the  ice-field  takes  place. 

The  one  numerical  comparison  made  shows  surprisingly  good  agreement 
of  theory  with  observation.  It  would  seem  that  a  study  of  more  realistic 
thickness  distributions  would  prove  fruitful. 


21 


BLANK  PAGE 


R- 13 13 


PART  II 

A.  MODEL  D.  WAVE  TRANSMISSION  THROUGH  A  FLEXIBLE  ICE-FIELD 

I.  FORMULATION 


SKETCH  2 


A  semi-infinite  ice  sheet  is  floating  on  the  surface  of  water  of 
constant  depth  H  .  A  plane  wave  is  obliquely  incident  from  the  region 
_  <»  <  y  <  o  .  Let  $(x,y,z,t)  he  the  velocity  potential  of  the  liquid 
motion.  Then  on  the  linearized  theory  of  smali  amplitude  water  waves,  $ 
satisfies 


§  +  §  +  §  =0 

0  <  z  <  H 

(LI) 

xx  yy  zz 

§  =  0 
z 

,  z  =  0 

(1.2) 

5  +  q®  -  0 

tt  9  z 

at  z  =  H 

(1.3) 

«,  CD 

<  y  <  0  (free  surface) 

Let  the  equation  of  the  ice  field  be 


z  =  §(x,y,t)  +  H 


(1.4) 
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where  l(x,y,t)  represents  the  displacements  of  the  ice  sheet  above  the 
undisturbed  free  surface  z  =  H  .  Then  Bernoulli's  equation  gives 

P  -  PQ  =  "  9?)  0-5) 

and  we  also  have 

as  _  af 

Sz  dt 

on  z  =  H  ,  on  the  linear  theory. 

It  Is  assumed  that  the  ice  sheet,  having  constant  mass  thickness  h  , 
is  displaced  from  equilibrium  by  the  differential  pressure  p  -  p  ,  and 
that  each  element  will  be  subjected  to  a  force  arising  from  the  bending 
stresses  in  the  sheet. 

It  may  be  shown ^  that 


(1-6) 


p  -  pQ  =  Dii4  §  +  p.h(Stt  +  g) 


0.7) 


where 


and 


D 


Eh3 

12(l-vs) 


E  =  Young's  modulus 
v  =  Poisson  1 s  rat io 
p.  =  density  of  ice 


S7* 


+ 


dx3 


Combining  (1.5)  and  (1.7)  and  differentiating  with  respect  to  t 


gives 


n  k  •" 

I  +  a$  =  _  —  v4®  -  — —  9 

tt  3  z  p  z  p  ttz 


z  =  H  ,  0  <  y  <  ®  (1.8) 
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where  (1.6)  has  been  used. 

Since  the  ice  extends  to  infinity  in  either  x-direction,  it  is 
possible  to  Subtract  out  the  x-dependence.  Assuming  also  a  time  harmonic 
dependence  of  frequency  u)  ,  let 

£ (x ,y ,2 , t)  =  Re  |  0(y,z)  ellcx  e  ,aJtj  (1-9) 

whence  0(y,z)  satisfies 


-2-®  +  —  -  k30  =  0  ,  0<z<H  ,  -  =>  <  y  <  ®  (1.10) 

3ys  a2* 

■f®  =  0  ,  z  =  0  ,  -®<y<®  (1.11) 


K0 


_  CO  <  y  <  0 


o  <  y  <  » 


where  K  =  w3/g  =  — r-  ,  M  =  D/pg 


.Kh 

j _ 

P 


(1.12) 


(1-13) 


Conditions  (1.10)  to  (1.13)  are  not  by  themselves  sufficient  to  give 
a  unique  solution  0(y,z)  .  Additional  conditions  regarding  the  vanishing 
of  the  bending  moment  and  shearing  force  at  the  edge  of  the  ice  field  will 
be  imposed  together  with  assumptions  regarding  the  form  of  the  solution 
for  y  =  .  Assumptions  concerning  the  behavior  of  0  and  its  deriva¬ 

tives  near  y  =  0  ,  z  =  H  which  ensure  that  Fourier  transforms  converge 
will  be  made  during  the  course  of  the  analysis.  These  assumptions  may  be 
verified  once  the  final  solution  is  obtained. 

Note  that  the  particular  case  k  =  0  corresponds  to  a  plane  wave 
which  is  normally  incident  upon  the  ice  field;  that  is,  whose  crests  are 
para  llel  to  y  =  0  . 
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2.  PRELIMINARY  DISCUSSION  OF  THE  SOLUTION 

The  eigenfunctions  of  Eqs.  (i.10),  (l.Il)  and  (1.12)  for  y  <  0  are 

exp  [±  y  (kS  +  a^)2]  cos  anz(n  =  1,2,  ...)  where  the  an  are  the 
roots  of  the  equation, 

K  cos  aH+asinaH-0  (2.1) 

n  n  n 

It  is  shown  in  Appendix  A  that  there  are  an  infinite  number  of  real 

roots  a  such  that  ja  |  <  la  .  |  ,  (n  =  *1,  ±2,  ...)  and  that  for 

n  n  *  n+i 

K>  0  there  are  also  two  pure  imaginary  roots  ±  ia  with  eigenfunctions 

x 

exp  [iy  ( k3  -  ao  )2]  cosh  aQz  .  There  are  no  other  roots.  For  a  plane 
wave,  k  <  aQ  so  that  there  may  be  propagation  in  either  y-direction  for 
y  <  0  ,  since  the  exponent  is  purely  imaginary.  Thus  for  y  negative, 
the  expected  form  of  a  bounded  solution  would  be 

X 

exp  [±  iy  (aoS  -  k3)"]  cosh  aQZ  +  0  (exp  ky) 


For  y  >  0  ,  the  situation  is  more  complicated.  The  eigenfunctions  of 
Eqs.  (1.10),  ( 1 . 1 1 )  and  (1.13)  are  of  the  form 


exp  [±y  (k2  +  bn  )  3  cos  b^2  (n  =  1,2,  ...) 


where  the  b  are  the  roots  of  the  eauation 
n 


K  cos  b  H  +  b  (l-L+Mb4)sinbH  =0  (2.2) 

n  n  \  n  /  n 


A  detailed  examination  of  the  roots  of  this  equation  is  made  in  Appendix  A. 

It  is  found  that  for  L  <  1  ,  M  >  0  ,  Eq.  (2.2)  has  two  purely  imagi¬ 

nary  roots  ±  ;bQ  ,  a  doubly  infinite  sequence  of  real  roots  *  t>n(n  =  1,2,  .. 
and  four  complex  roots  ±  (the  assumption  L  <  1  covers  the  range 

of  practical  interest;  the  case  L>  1  is  not  considered  here).  The  exact 
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position  of  the  complex  roots  is  not  required.  It  is  sufficient  for  the 

time  being  to  note  that  c  =  \  +  ip.  where  A.  >  p,  >  0  ;  this  is  proved 

o 

in  Appendix  A.  Because  of  this,  c  2  lies  in  the  first  quadrant,  and  ^ 

Re  c  3  >  0  .  Then  it  follows  that  Re  c  1  >  k  where  c  1  =  (c  3  +  k3)2 
o  o  o  o 

the  square  root  being  such  that  c  '  =  c^  when  k  =  0  .  A  similar  ine- 

x 

quality  holds  for  c  1  =  (c  3  +  k3)2  .  Thus  for  y  >  0  there  exist 
o  o 

eigenfunctions  of  the  form 

exp  |  ±y  (k3  -  b  V  cosh  b  z 


2  2 

which  provide  the  wave  propagation  if  >  k  , 


exp ±y  (k3  +  bn3)~|  cosh  b^z  ,  (n  =  1,2,  — )  , 


and 


exp  f ±y  c^1  >  cosh  c^z 


exp  •{  iy  cQ'  t  cosh  cqz 


Thus  for  y  positive,  the  expected  form  of  a  bounded  solution  would  be 


expj^iy  (b^3  -  k3)2  j  cosh  bQz  +  0(e 

if  b  2  >  k3  ,  the  terms  0(e  ky)  arising  from  the  fact  that  Re  c  '  = 
o  o 

Re  c  •  >  k  ,  If  b^3  <  ks  ,  then  there  will  be  no  wave  propagation  into 

the  ice.  This  is  made  clear  in  Fig.  2.  For  the  incident  wave,  we  write 

k  =  a  sin  9  ,  and  for  the  transmitted  wave,  k  =  b  sin  9*.  where  6_ 
o’  o  T  T 

is  the  angle  of  transmission. 

Then  the  incident  wave  is  of  the  form 


i a  y  cos  S  +  ia  x  sin  9 
o  o  , 

e  cosh  a  z 

o 


whereas  the  transmitted  wave  takes  the  form 
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Clearly,  since 


Eb  y 
o 


0.  +  ib  x  sin  9t 

To  T 


cosh 


b  z 
o 


(2.3) 


propagation  into  the  ice  will  only  occur  if  sin  9  <  .  Whenever 

aQ  sin  9  >  ,  the  transmitted  wave  becomes  exponentially  damped  of  the 

form 

1 

-y(a  2  sin2  9-b  2) 

e  °  0  cosh  b  2 

o 


so  that  in  this  case  a  1 1  modes  decay  exponentially  and  the  incident  wave 
is  totally  reflected.  Thus  for  given  K,  L,  M,  there  exists  a  critical 
incident  angle  ecrjt  such  that  an  incident  wave  approaching  the  ice-field 
at  an  angle  9  >  9cr^  is  completely  reflected  (see  Fig.  3)«  This  will 
occur  when  9y  =  S0°  so  that 


0  . 
cr  1 1 


sin  (b  /a  ) 
o  o 


(2.4) 


The  transmitted  angle  9^.  indicates  whether  a  given  incident  wave 
will  be  bent  towards  or  away  from  the  normal  to  the  ice  field.  If 
b^  <  then  from  Eg.  (2.3),  @  <  9^.  ,  the  transmitted  wave  is  bent  away 

from  the  normal,  and  there  will  always  exist  a  critical  angle  given  by 
Eq.  (2.4).  On  the  other  hand,  if  a^  <  ,  then  6^.  <  6  and  the  trans¬ 

mitted  wave  is  bent  towards  the  normal  (see  Fig.  4).  In  this  case,  an 
incident  wave  approaching  the  ice-field  will  always  penetrate  the  ice 

regardless  of  the  incident  angle.  If  a  -  b  ,  no  deviation  of  the 

o  o 

incident  wave  occurs. 

For  the  time  beinq,  it  will  be  assumed  that  a  sin  6  <  b  so  that 

o  o 

there  exists  an  undamped  progressive  wave  travelling  into  the  ice-field. 
In  the  case  of  normal  incidence,  9=0,  there  will  always  exist  such  a 
wave. 
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3.  METHOD  OF  SOLUTION 


(a)  Derivation  of  the  Wiener-Hopf  Equation 

The  solution  for  the  function  0(y,z)  is  achieved  by  means  of 
Fourier  transforms  and  the  Wiener-Hopf  technique.  In  order  that  the 
Fourier  transforms  might  converge  in  a  strip  of  the  transform  variable 
plane,  the  following  device  is  used.  The  preceding  section  indicates  the 
expected  form  of  the  solution  for  large  values  of  y  .  It  is  anticipated 
that  a  prescribed  oblique  plane  wave  incident  from  y  <  0  will  give  rise 
to  a  reflected  wave  in  y  <  0  ,  and  a  transmitted  wave  in  y  >  0  .  The 
amplitude  and  phase  of  the  reflected  and  transmitted  waves  are  determined 
once  the  incident  wave  is  prescribed.  However,  we  shall  fix  the  amplitude 
and  phase  of  the  transmitted  wave  beforehand  and,  hence,  determine  the 
reflected  and  incident  waves.  The  reason  for  this  will  soon  become  apparent. 


Thus,  let 


i  fa  'y 

0(y,z)  =  e  0  cosh  bQZ  +  0(e  y)  ,  O^z^H  ,y>0 


(3.0 


where  b  '  =  (b  8  -  k2)2  and  b  1  =  b  when  k  =  0 
o  o  o  o 


■ia  y  -ia  'y 

and  0(y,z)  =  Ie  0  cosh  a  z  +  Re  °  cosh  a  z  +  0(e 


O^zSH  ,  y  <  0 

(3.2) 

where 

a  1  =  (a  2  -  k2)~  and  a  1  =  a  when  k  =  0 
o  o  o  o 


Division  of  the  solution  by  I  gives  the  solution  due  to  an  inci- 
ia0[y 

dent  potential  e  cosh  a  z  . 

o 

Consider,  now,  the  function  ’Ky.z)  where 


ib  !y 

t(y,z)  =  0(y,z)  -  e  °  cosh  b  z 

o 
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Then  from  (3.1)  Ky>z)  is  exponentially  small  for  y  >  0  so  that 
the  Fourier  transform  of  V(y,z)  with  respect  to  y  will  exist  in  a 
strip  of  the  transform  plane;  a  basic  requirement  for  the  successful  appli¬ 
cation  of  the  Wiener-Hopf  technique. 

Now,  t(y,z)  satisfies 

^4  +  ”-tc%  =  0p0<z<H  ,  -»  <  y  <  «  (3.3) 

£>y3  a2 

|i=0,z«=0,-»<y<«  (3.4) 

.  ib'y 

KV=-r^-  +  Ae  °  ,z  =  H,-«=<y<0 


where  A  =  (Mb0*  -  L)bQsinh  bQH 


Kt  =  (1-L)  |i+  M(^4  -  k2)  IJ  .  z  =  H,0<y<-  (3.6) 

6y 

My.z)  =  0(e  ky)  ,  y  >0  ,  for  each  z  (3.7) 

+  Ia'y  -ia*y  ib‘y  , 

+  (y,z)  =  Ie  cosh  a  z  +  Re  °  cosh  a  z-e  °  cosh  b  z  +  0(e  , 

o  o  o 

y  <  0  for  each  z  (3.8) 

It  is  assumed  that  outside  some  neighborhood  of  (0,H)  (i  and  its 
first  and  second  partial  derivatives  are  also 

0(e"ky)  ,  y>  0  and  0(1)  ,  y  <0  (3-9) 
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It  is  further  assumed  that  if  is  bounded  everywhere  and  that  in  a  neigh¬ 
borhood  of  (0,K) 


5i|i  5ij' 

by  ’  dz 


0  <  0  <  I 


(3-10) 


where  r2  =  y2  +  (z  -  H)a 

The  reason  for  the  assumption  given  by  Eq.  (3.10)  requires  some 
explanation.  There  is  no  a  priori  reason  why  the  velocity  components 
should  be  non-singular  at  the  edge  of  the  ice-field.  Such  s ingu lar it ies 
invariably  occur  in  potential  problems  at  the  confluence  of  two  boundaries 
on  which  different  boundary  conditions  are  satisfied.  However,  physically, 
we  require  that  there  be  no  breaking  of  waves  at  the  interface  as  this 
would  introduce  an  arbitrary  constant  to  determine  the  amount  of  energy 
loss  which  occurs.  On  a  linear  theory  such  a  breaking  phenomenon  is 
represented  by  a  sink  singularity  in  i|r  so  that  the  velocity  components 
would  be  0(— )  in  the  neighborhood  of  the  leading  edge,  corresponding  to 
the  logarithmic  singularity  in  i|>  .  By  insisting  that  the  singularity  be 
of  the  form  (3.10),  loss  of  energy  due  to  breaking  is  just  avoided,  at 
the  same  time  allowing  any  milder  singularity  to  occur  in  the  analysis. 

From  conditions  (3.9)  and  (3.i0),  the  Fourier  transform 


*(or,z) 


(3.11) 


exists  for  0  s  z  ^  H  ,  and  is  regular  for  o>  (=  a  +  It)  ;n  the  strip 
D:  -  k  <  t  <  0  of  the  complex  <*-plane  (see  Sketch  3,  p.  37).  From 
Eq-  (3.3) 

,a3'f  5a  V  .3,,  iay,  n 
(~r  +  —  -  k  $)e  'dy  =  0 

3y  9z3 

and  integration  by  parts  gives 

["  -  (a2  +  k2)  1  f(o',z)  =  C  ,  (cr  e  D  ,  0  £  z  <  H) 

L  a, 3  -1 
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Define  Y  =  (o3  +  k2)2  such  that  the  cuts  extend  from  ±  ik  to 
along  the  imaginary  axis  in  the  a-plane.  Then, 

(i)  As  «  =  c  (real)  -  »  ,  Y(cr)  -  c 

(ii)  Re  (y)  >  0  everywhere  away  from  the  cuts 

and 

T(a,z)  =  Ax (a)  cosh  yz  +  A2(a)  sinh  Yz 
The  transform  of  condition  (3.^)  is 


1  im 
z—  0 


ST ,  . 


o 


so  that  As  (a)  =  0  . 

From  condition  (3.10)  it  may  be  shown  that 


1  im 
z  -'H  - 


and 

1  im 
z  -*H- 


Thus 


S  <*•*>  =  I 


and  (a,H)  =  YH  tanh  YH  V  (a.K)  (»  £  D) 

In  order  to  obtain  a  Wiener-Hopf  equation,  the  following  notation 
Let 

Hat. z)  *  »+(<»,*)  +  T.fcf.z) 


(3.12) 

(3.13) 

is  used. 
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where 


exists  for  0  ^  z  ^  H 
and 


r+(cr.z)«J  Ky.z)e!ffy 


dy 


and  is  a  regular  function  of  o 


in  0+: 


'Ky.z)e‘£Vydy 


exists  for  0  s:  z  ^  H  and  is  a  regular  function  of  or  in  D  : 
Then 


lim  Y,  {a,z)  =  Y  (o,h) 


z-H- 


3Y± 


3Y., 


lim  ~dz  =  "ST"  (a»H) 

z  H  - 


3Ya 


In  future  Y±  (c*,H)  and  (a,H)  will  be  abbreviated  to 

res pectively. 

Now  the  transform  of  condition  (3-5)  is  easily  seen  to  be 


KY_  =  Yf. 


iA 

w  +  bo 


(a  e  D_) 


Condition  (3-5)  is  more  troublesome.  It  is  assumed  that  all  y 
dii  ,  .  ^  . 

of  (y,H)  up  to  and  including  y,-4-g~z  (y,H)  are 


0(e"ky)  ,  y  >  0 


It  is  also  assumed  that 


(y,H)  is  o <-L) 
dy  Sz  y 


(0  <  P  <  l )  as  y  -  0+ 


(3.14) 
>  -  k 


(3-15) 

<  0 


and  Y±‘  , 


(3.16) 

derivatives 


(3.17) 


(3.18) 
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Then  the  transform  of  condition  (3.6)  exists  and  integration  by  parts 
gives 


%so 

'4> 


dy  =  -  04  +  i^cs  +  +  2k'J)(cs  -  iaey) 


where  condition  (3-17)  has  been  used,  and  where 


+  y4  Jf+>  (a  s  D+) 

(3.19) 


y=0+ 


a 

Cz~W 


[M. 


-  0+ 


6y  Jy=0+  ^  dy° 

Thus  the  transform  of  condition  (3.16)  becomes 


y=0+ 


KY+  =  ¥  '  [1  -  L  +  MY4]  -  M  [c^-ioe,  -  (o2  +  2ka)  (o,  -  i  a  ct )  ]  (a  s  D+) 

(3.20) 

Now  add  Eq.  (3.16)  to  Eg.  (3 - 20)  and  use  Eq.  (3.13).  Then 
K  cosh  yH  -  Y  ( 1  -  L  +  My4)  sinh  yH  j  Y+'  +  j  K  cosh  YH  -  Y  sinh  YH 

=  -  Y  sinh  YH  j  -^-v  +  M  [c4-i»c3  -  (o?  +  2k2)  (Cg-ictc,)]  J  (a  e  D) 

(3.21) 

Equation  (3.21)  is  a  typical  Wiener-Hopf  equation  holding  in  a  strip  of 
the  complex  «-plane.  Although  Eq.  (3.21)  involves  the  constants  c.(i=T,2, 
3,^),  it  will  be  shown  how  f'  and  hence  T  may  be  determined  in  terms 
of  just  two  constants. 

Let 

K  cosh  YH  -  y(1  -  L  +  MY4)  sinh  YH  =  f  (a) 

(3.22) 

K  cosh  yH  -  Y  sinh  YH  =  f  (a) 

o 
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Then 


f  -  f  =  -  (MY*  -  L)  Y  sinh  YH 
1  o 


(3-23) 


and  Eq.  (3-21)  becomes 


fx  (My4  -  L)  ¥  ■  -  ,  -  M  rc,-io-c3  -  (<^+  2k2)(ca-;acJ] 

o  J 

-  -  fDJ(MY4  -  L)  »_■  +  -  |Ab  y  +  M  [c4-i^c3  -  (°?  +  2k2)(c2-iaCl]J  (ore  0) 


(3.24) 


It  is  shown  in  Appendix  B  that  we  may  define 


f1(»)  _  K+(a) 

f  (a)  ~  K  (a) 
o  — 


(3-25) 


where  K±(a)  is  regular  and  non-zero  in  D±  ,  respectively.  Furthermore, 
it  is  shown  that 


and 


+* 

f — \ 

£ 

II 

o 

s 

% 

V 

|ce|  —  “  ,  (a  c  D^_) 

(3.26) 

K_(“)  =  0(-^)  , 

i 

o 

(j 

8 

t 

~ 

(3.27) 

Then  (3.24)  may  be  written 


(  1  fK  (o')  -  K  (-  b  1  ] 

K+(»)  (My'  -  L)  *+‘  -  M  tc4-iac3  -  (a2  +  2k2)  (c^-IOc,  )]  -  iA  — - —  ,  --p 

*  *  L  o 


*  -  K_(a)  J  (My4  -  L)  Y_*  +  M  [c^-iarc3 


(ti2  +2  k2  )  ( c2  -  i  OCj  ] 


I 

i 


-  iA 


'*_(“)  -  *+<"  b0’)' 
a  +  b  ' 

*-  o-i 


(a  e  D) 

(3.28) 
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It  is  clear  that  each  term  on  the  left-hand  side  of  (3-28)  is  regular  for 
<*  e  D+  ,  whereas  each  term  on  the  right-hand  side  is  regular  for  or  €  D 
Since  the  two  sides  are  equal  for  a  e  D  ,  this  defines  a  function  j(o) 
regular  for  all  (finite)  a  .  The  determination  of  J(a)  is  achieved  by 
a  consideration  of  its  behavior  as  |aj  —  <*>  . 

Consider  the  left-hand  side  of  Eg.  (3-28).  From  (3.20),  it  is  equal 
to 


K+wfe-n}  - 


fK+(a)  -M-b^) 

iftl  ^ 


(a  e  D+) 


Now  ¥  ,  —  0  as  |a|  “  m  in  D+  from  condition  (3-10),  (see  Noble, ^ 

p.  38,  Eq.  fl.74])." 


Thus,  since  K  (a)  =  0(0^)  ,  ja|  -*  ®  in  D+  ,  then  J(a)  =  o(as)  as 

*  -I 

| q,|  _  oo  jn  p+  .  Similarly,  since  K_(or)  =  OOjnr)  as  j rv|  — »  00  in  D_  , 

and  1  '  —  0  as  |<*|  in  D  ,  then  J(a)  =  ot0^)  as  | cr |  -»  ®  in  0_  . 


Hence 


J («)  =  0(0^)  as  \a\  -  ® 


(3.29) 


Thus  by  an  extension  of  Liouvi lie's  theorem, 


(Sec.  2.52). 


J(ff)  =  CjO  +  Cg  (CL ,  C3  constants)  (3.30) 


Equating  each  side  of  Eq .  (3.28)  to  JO*)  ,  solving  for  ¥  1  and  Y 
and  adding,  gives 


(h-y4-l)  +  *:) 


(M^-L)Y'  -  (j(a) 


iAK^-b;Hf  1  1  , 

Oi+b^  J  “  iGSfr 


Note  that  this  step  is  crucial  to  the  ensuing  argument.  If  we  let  3  =  1 
In  condition  (3-10)  indicating  an  energy  source  at  the  origin,  then  all  we 
can  say  is  that  ,  Y  1  are  bounded  as  |o|  —  “  in  0+  and  in  the  sub¬ 
sequent  argument  we  find  that  J(aO  =  +  C,_®  +  C3  .  Then  the  solution 

is  no  longer  unique;  the  strength  of  the  energy  source  at  the  leading  edge 
of  the  ice-field  must  be  given  i r>  order  to  determine  the  additional  constant. 
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From  Eq.  (3.13),  this  may  be  written  in  terms  of  ¥  .  Thus, 


r.  !AK+(-bO,  . 

v(c.h)  -  - -1 - r 

°  K+(a)|K-y  tsnh  yHj- 


(3.31) 


V  (a,H)  ■=>  {j^)  -  — 


K.M^K-yO-L+My^)  tanh  yHj- 


(3-32) 


Equation  (3.12)  and  the  Fourier  inversion  formula  give 


y/v  z)  =  _L  f  Y(a  H)  -0sh  -~V2  e"'aydw 

V”‘  J  2TT  Jc  cosh  yH  e  da 


(3.33) 


where  C  is  some  path  in  D  as  shown  in  Sketch  3- 


SKETCH  3 


a-PLANE 


-a0  -b0 


.  '  1 
ao 
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The  constants  Ct  and  C2  will  be  determined  from  the  conditions  of  zero 
shearing  force  and  bending  moment  at  the  edge  of  the  ice,  z  =  H  ,  y  =  0+  . 
Assume  for  the  time  being  that  this  has  been  done.  Then  it  is  necessary  to 
check  that  Eq.  (3.33)  with  ’*'(<*, H)  ,  given  by  either  Eq.  (3-31)  or  Eq .  (3-32) 
does  in  fact  satisfy  the  conditions  of  the  problem.  Now  T(ar,h)  is  0(*'”s) 
as  |<z|  -*  m  in  D  so  that  the  integral  in  Eq.  (3.33)  is  uniformly  con¬ 
vergent  for  0  £  z  ^  H  ,  all  y.  Also,  for  0  s  z  <  H  ,  all  y  ,  the 
integrals  obtained  by  differentiating  |(y,2)  with  respect  to  y  or  z  , 
any  number  of  times  are  also  uniformly  convergent.  Since  the  operator 

-  k2  applied  to  f(y,z)  makes  the  integrand  vanish  identically, 

it  is  clear  that  condition  (3.3)  is  satisfied  for  0  £  z  <  H  ,  all  y  . 

Similarly,  condition  (3.A-)  is  seen  to  be  satisfied  on  z  =  0  .  To  verify 

condition  (3.5)  for  y  <  0  ,  Eq.  (3.31)  is  used  in  Eq .  (3.33).  and  the  path 

of  integration  is  deformed  upwards  in  such  a  way  that  the  ends  of  the  path 

tend  to  infinity  along  lines  in  the  upper  half-plane.  It  is  now  permissible 
to  apply  the  operator  K  -  to  iKy,z)  and  put  z  =  H  .  Thus 


(K* 


•AK+(-b^)  ^  e-idy 

B+b'  J  K+(a) 
o 


da 


where  C’  is  the  deformed  path.  It  is  clear  that  the  only  contribution 
to  the  integral  arises  from  the  only  singularity  of  the  integrand  in  D+  , 
which  is  a  pole  at  a  -  -  b  1  .  Thus,  for  y  <  0  , 


W  -§£) 

02  z=H 


Ae 


iboy 


which  verifies  condition  (3.5). 

To  verify  condition  (3.6)  for  y  >  0  ,  Eq.  (3.32)  is  used  in  Eq.  (3.33), 
and  the  ends  of  the  path  of  integration  are  deformed  downwards  so  as  to 
tend  to  infinity  along  lines  outside  D  in  the  lower  half-plane.  Now  for 
y  >  0  ,  it  is  possible  to  differentiate  |(y,z)  any  number  of  times  and 
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put  2  =  H  .  Thus 


{k  -  0-l)  -  '<*)  'Ky.H) 

dy 


iAK+C-fa*)  .  -<ay 

“W^-/fcRda 


=  o 


The  boundedness  of  t(y,z)  everywhere  follows  from  the  fact  that  +(y,z) 
is  given  by  a  uniformly  convergent  integral  for  0  s  z  ^  H  ,  all  y.  A 
differentiation  with  respect  to  either  y  or  z  produces  an  additional  a 
into  the  numerator  of  the  integrand  so  that  the  discontinuous  part  of  the 
integral  behaves  in  a  similar  manner  to 


e-iorye-w(H-z) 

i  =  J  - -  —  -j -  dct  ~  log  r  +  bounded  terms  near  z  =  H  ,  y  =  0 

Thus,  condition  (3-10)  is  verified. 

Condition  (3-7)  is  verified  by  using  Eq.  (3.32)  in  Eq.  (3.33)  and 
deforming  the  path  of  integration  into  the  lower  half-plans,  taking  into 
account  the  poles  of  the  integrand.  Since  the  singularity  of  the  integrand 
nearest  the  real  ff-axis  is  at  a  distance  greater  than  k  from  it,  it  is 
clear  that  ilr  and  also  its  first  and  second  partial  derivatives  are 
0(e  ^7)  >  f0r  y  >  0  and  each  z  .  A  similar  argument  may  be  used  to 
verify  condition  (3.8).  In  this  case,  contributions  are  obtained  from 
the  poles  of  the  integrand  at  a  =  -  bQ'  ,  a  -  ±  a^1  ,  together  with 
contributions  of  0(e  ^)  . 

The  conditions  (3.17)  and  (3-l8)  require  special  attention.  Now 


<_(<*)  [k-yO- 


•igyYtanh  YH 


L+My4) tanh 


dor 


YH 


3S 
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where 


Fj  (a)  =  J  (a)  - 


!AK+(-b;) 

**o 


and  where  the  path  C  has  been  deformed  so  as  to  pass  beneath  the  point 
a  =  -  ik  .  Write 


Y  tan  h  Mi 


1 


K-y(l-L)tanh  ~yti 


K-Y(l-L+My4)tanh  yH  My*  L  K-yO -L+My*  )tanh  yH 


-•} 


Then 


\  («) 

HY4K  (a) 


is  regular  and  tends  to  zero 


Im  (o)  <  -  k  ,  so  that 


j-  r 

2tt  4 


Fl  (a)  e 
My4  K_(a) 


-lay 


dc  =  0 


Thus 


1  p  Fj  (a)e  '“y-[K-Y(l-L)tanh  yij- 
2TT  Hy4 K_  («){k-y ( 1  -L+fly4 )  tanh  Yh} 


(3-34) 


where  the  integrand  is 


0<£> 


It  is  clear  that  all  y  derivatives  up  to  and  including  — nr(T^) 

°y  dz  H 

are  bounded  for  y  ^  0  ,  and  deforming  the  path  of  integration  around 
the  poles  of  the  integrand  indicates  that  all  derivatives  are  0(e-l<y) 
for  y  >  0  .  It  may  be  shown  in  a  similar  manner  to  the  verification  of 
condition  (3.1)  that 


H 


O(log  r) 


r2  =  y2  +  (z-H)2  near  r  =  0 


Thus  conditions  (3-17)  and  (3.18)  are  verified,  and  t(y,z)  does  indeed 
satisfy  the  conditions  of  the  problem. 
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(b)  The  Far  Field 

In  order  to  demonstrate  the  form  of  the  solution  for  large  y  ,  it 
is  convenient  to  use  series  representation  of  the  solution. 

Thus,  for  y  >  0  ,  the  integral  expression  (3.33)  may  be  written 
in  terms  of  a  sum  of  the  residues  at  the  poles  of  the  integrands  (see 
Appendix  B) .  Thus, 


Ky.z) 


-c;y 

2F,{-ic,)e  c  cose  z  cose  H 
_ 1  o _ o  o _ o _ 

co[2 co  H0 -L+Mco4)+(  1  -L^5Mc04 ) s in2 cq h]k_ (-id) 


2F,  (-1  c1  )eC°y c  cosc  z  cose  K 

1  o _ o  o _ o _ 

c^nO-LttTc;  )+(1-L+5MCq4  ) s i n 2 h|k_ (- i  e^) 

-biv 


2Fi  (-ib^)e 


b  cos  b  2  cos  b  H 
n  n  n 


^7  p  1 

n=T  b‘  2b  H (l-L+Hb4)+(l-L+5Hb4)sin2b  H  [K  (-lb1) 
nL  n  n  n  n  J  ~  n 


For  y  <  0  ,  similarly, 
ib  'y 

>My,z)  =  -  e  °  cosh  bQz 


(3.35) 


2ia  cosh  a  z  cosh  a  H 
a  a  o 


a1  (2s  H  +  si  nh  2a  H) 
oo  o 


CO 

-E 

n=l 


-1a1 y  ia‘y 

f  Fi(ao)e  °  Fl(~ao^e  0 


K+  <*;>> 


K+ 


b'y 

2F.  (ibnOe  b  cosb  z  cos  b  H 
_ n  _ n  n 


b  '  |2b  H  +  sin  2b  h] 
n  l  n  n  J 
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Comparison  of  the  coefficients  In  Eqs.  (3-35)  and  (3.36)  with  the 
coefficients  ft  and  I  occurring  in  Eqs.  (3.1)  and  (3.2)  indicates  that 


and 


r  M+  V>  V-  V> 

I  '  M"  K+<+  aQ') 


(3.37) 


(X)"1 


i a  '  (2a  H  +  sin  h  2a  H) 

O  O  O 

2a  cosh  a  H 
o  o 


M-  V> 


(3.38) 


The  reflection  coefficient  &  is  defined  as  the  ratio  of  the  amplitude  of 
the  reflected  wave  to  the  incident  wave  at  infinity.  The  transmission 
coefficient  "T  is  defined  similarly. 

Now 


d£  _ 

dt  ^dz' 


z=H 


so  that  the  elevation  §  is  given  by 


l(x.y.t)  = 


i  kx-iiutl 
e  ; 


Thus  the  elevation  of  the  incident  wave  is 


C(x,y,t) 


fa  sinh  ®  H  ikx+  iaj y  -  iu>ty 

»»{— -i.—  1 «  °  } 


(3.39) 


with  ampl itude 


I 


a  sinh 
o 


a  H 
o 


(3.40) 


42 


R-1313 


Similarly  the  amplitude  of  the  reflected  wave  is 


a  sinh  a  H 

|r|  -g— . -°  ■  (3.41) 

1  1  (D 


The  elevation  of  the  transmitted  wave  is 


!b  sinh  b  H  i kx  +  i b  *y -  i 1 
~ - —  e  °  |  (3.42) 

b  sinh  b  H 

with  amplitude  — - —  so  that  the  reflection  coefficient  is 


a 


1 R 1  _ 
II! 


W=^T\ 


(3.43) 


since 


and  the  transmission  coefficient  is 


(see  Appendix  B) 


r 


b  sinh  b  H 

o _ o_ 

a  sinh  a  H 
o  o 


III 


^(2V+sinh  2aoH)s?nh  boH  |K+(-a;)| 

a  3  sinh  2 a  H  !  ( “ao^  I 

o  o 

(3.44) 


(c)  Determination  of  J(a) 

There  are  still  two  conditions  which  need  to  be  satisfied  at  the 

leading  edge  of  the  ice  before  the  solution  is  uniquely  determined.  These 

are  conditions  which  express  the  fact  that  no  energy  is  put  in  or  taken 

out  of  the  system  at  the  leading  edge  of  the  ice-field.  A  derivation  of 

13 

these  conditions  is  given  by  Hildebrand.  They  are 
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^2  -  2 

Oy  dx 


(3.45) 


oL(^  +  ^l)  +  0O-v)-A 


dy  3yS  dx2 


(3.46) 


BySx* 


where 


_  a  /’a®?  sEA  . 

D  3y  [dy*  +  dx3  j  'S  the 


(1-W  D 


Sag 

dydx 


transverse  shearing  force  along  y  =  0 

is  the  twisting  moment  along  y  =  0 

is  the  bending  moment  along  y  =  0 


We  have 


i  =  Re 


9 


o0  ikx-iu)t 


i  in  3z 


and  the  conditions  (3-45)  and  (3-46)  in  terms  of  I  ,  become 


o2  ,3j5 


—G*)  -  vk~  ££}  =0  ,  y  =  0+ 

a.. a  °Z  Z=H  OZ/ 


ay* 


z  z=H 


(3.47) 


a3  /Bjzi 


s— (SS.)  (2-v)  k3  |-(|^)  =  0  ,  y  -=.0+ 

5y3voz'z"H  oy  ',dz'z=H  1 


(3.48) 


These  conditions,  in  turn,  may  be  written  in  terms  of 


_H-(2-v)k®  „  =  ib  b'sinh  b  H 

^y3  oz'z=H  '  oy'oz',z=H  o  o  o 


=  ib  b'sinh  b_H 
o  o  o 


W(y,z).  Thus 

b^+(2-v)k2 
f>o®  +  (l-v)k2 


(3-49) 
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and 


5s 


by5 


02=H  ~  Vk2(^z=H  =  bo  iinh  b«H  [bo'=  +  vk“] 


=  b^  slnh  bQH  jb^2  -  ()-v)k2J  y  =  Of 

(3-50) 

In  the  case  of  normal  incidence,  k  =  0  ,  and  the  conditions  on  0  are  just 


~(|f)  =  0  ,  y  =  Of 

3y  32  z=H 


(3-51) 


d3  ,30, 


“(£>  =  0  ,  y  =  0+ 

*y3  az  z=H 


(3.52) 


which  compare  with  Stoker 

Now  if  the  expression  (3-35)  for  t(y,z)  is  differentiated  with 
respect  to  z  we  obtain,  for  z  =  H  , 


kSz'z=H 


Fi  (-'co)e"C°y9(c0)+F1  (-'CQ)e-coyg(Eo) 


"2  Ft('5bn)e  n  9(bn) 

n=l 


where 


(3-53) 


g(e> 


2(K  cos  PH  f  g(l-L)  sin  PH)  cos  0H 
HP3  •  p'K  {-iP')  [2PH(  1-L+MP4)  +  (1-L+5MP4)  sin  2PH] 


and 

JL 

e*  =  +  (p3  -  k2) 2 


^5 
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Remember  that 


Fx  (a)  =  Cxa  +  C2  - 


iAK  (-  b  ■) 

+ _ Q 

(a+b  0 


The  conditions  (3.4-9)  and  (3.50)  when  applied  to  (3.53)  give 


(c^3 -vk3 ) F1  {-i c^)g  (cq)  +  -vk2  ) Fx  <-i c^) g  (cj 


o’"'  o' 


cn  _ 

+  ^  (b^2-vk3)F1  (-Ibn)g(bn)  =bosinh  b^H  -(l-v)k3J  (3-54) 

n=»  I 


as 

£  t>n[bn2-(2-V^2],;i(“ibn)3(t'n)  "  boH  (V  " 

n=l 


Now 


(3.55) 


iAK  (-  b  ■) 
F:  («)  =  C2  +  Cxtx  -  - a~ 


ar+b  1 
o 


where 


A  =  (Mb  4  -  L)  b  sinh  b  H 
o  o  o 


Cl  early,  there  exist  two  equations  for  the  determination  of  Cx  and  C2  . 
After  some  algebra  they  may  be  written  in  the  form 


An  Ci  +  AiS  C3  =  Bx 


A_,  C.  +  A  C  =  B 
21  1  22  2  2 


(3.56) 
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where 


A,  -  -  ic  ,s  tc  ,£  -  (2-v)  k2J  g(c  )  -  ic  ,s  tc  ,2  -  (2-v)  k2]  g(c  ) 


O  Q 


o  o  o 


1  bn*3  [bn‘2  -  (2"V)  ^  s(b) 
n  n  n 

n=l 


(3.57) 


Alz  =  cn'  \cja  -  (2-v)  k3]  g(cQ)  +  cq'  [cQ'3  -  (2-v)  k3]  g(cn) 


o  o 


+  2  bn'  [bn,a  ‘  (2'V)  ^  9(bn) 

n=l 


(3.58) 


8i  =  -  ibo  bo’  Sinh  boH  lho  -  (‘-V)  k*]  -  AK+(-  bo'>  c'+°ib>  [co'a-(2-v)k3]  g(co 

(  o  o 

c  '  °°  b  1 

+  -  ° . ,  ;  Tc'3  -  (2-v)  k2]  g(c  )  +  ^  ,  ,  n ..  (  (b  ,s  -  (2-v)k3)  g(b  ) 

c  1  +  tb  1  o  o'  /  -j  b  1  +  ib  1  n  n 


o  o 


n=  I 


Asi  =  -  ic«'(c^'2  -  vk3)  g(0  -  (c^tS- vkS)  g(£0) 


-  1  E  b„'<b„'*  -  *“>  3<b„> 

n=l 


(3  .55) 


(3.60) 


-  (co=2  -  vk3)  g(co)  +  (So  =  2-vk2)  g(5n) 


+  2  (bn'3-vk2)  g(b  ) 
n=I 


(3.61) 
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B  =  b  s  in  h  b  H  fb  3  -  (1-v)k2} 


-  AK  (-  b  ') 
+  o 


(C>2-Vk3)  (E'3-vkz) 

c  ■+  ;b  ■'  9<Co)  +  ■  "g  |— ~rb  .  g(CQ) 


+ 


^  (b  ’2-vk2} 
jLj  (b  ’  +  ib  ')  g(bn* 

n  =  l  n  ° 


(3.62) 


Once  Cj  and  Cs  have  been  determined,  the  solution  of  the  problem  is 
complete  and,  in  particular,  the  expressions  (3.43)  and  (3-44)  for  the 
reflection  and  transmission  coefficients  may  be  computed. 

Unfortunately,  the  determination  of  the  constants  C  and  C?  above 
presents  enormous  computational  difficulties,  involving  as  it  does  the 
complex  roots  b^  of  the  transcendental  equation 


K  cos  b  H  +  b  (1-L  +  Mb  4)  sin  b  H  =  0 
n  n  n  n 

Thus,  whereas  the  solution  to  the  full  linear  problem  has  been  obtained, 
it  will  be  necessary  to  resort  to  the  shallow-water  approximation  to 
derive  the  solution  in  a  form  more  suited  to  computation.  It  is  possible, 
however,  to  derive  a  relation  between  H  and  which  will  be  done  in  the 
following  section. 


4.  RELATION  BETWEEN  H  AND  T" 

Although  the  explicit  determination  of  and  presents  enormous 
computational  difficulties,  it  is  possible  to  derive  a  simple  analytical 
relation  between  ~/Z.  and  ,  by  an  application  of  Green's  theorem.  We 
have  seen  that  the  two-dimensional  function  0(y,z)  satisfies  the  follow¬ 
ing  conditions  (Eqs.  1.10-1.13). 
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$.  +  (Lii-k20  =  o.O<z<H1-®<:y<®  (4.1) 

-  E  2 

Oy  oz 


<  y  <  a> 


(4.2) 


KfS 


<  y  <  0 


(4.3) 


K0  =  (1-L)|f  +  M^s-  -k2)  ||  ,  z  =  H,  0  <  y  <  «  (4.4) 

Oy 

Also  0  is  continuous  and  satisfies  the  conditions  (3.47)  and  (3.48), 
namely 


S2  ,dj5. 

57fe 


z=H 


I  2 

vk  {^5 


=  0 


z=H 


y  =  0+ 


(4.5) 


a"*  ,a£k 


(2-v)  k: 


a  ,o£i, 

°y  2=h 


=  0  y  =  0+ 


(4.6) 


Green's  theorem  applied  to  the  functions  0(y,z)  and  its  conjugate 
0(y,z)  may  be  written 


0v20)dydz  =  Jc(0  ||  -  i  ||)ds  (4.7) 


where  S  is  the  area  of  the  rectangle  formed  by  the  lines  z  =  0  ,  z  =  H  , 
and  y  =  ~'iQ  >  C  is  the  boundary  of  this  rectangle,  and  n  denotes  the 
normal  directed  to  the  exterior  of  the  boundary.  Since  0  and  hence  0 
satisfy  Eq.  (4.1),  the  left-hand  side  of  (4.7)  vanishes  and  the  equation 
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reduces  to 


Im  Jc  s  sk  ds  " 0 


(4.8) 


It  Is  assumed  that  yQ  is  so  large  that  for  y  =  ±yQ  >  0(y , z)  may  be 
replaced  by  its  asymptotic  form  as  given  by  Eqs. (3-1)  and  (3.2).  Thus 


ib  'y 

0(y  ,z)  ~  e  °  °  cosh  b  z  ,  Q  £  z  ^  H 

wo  c 


(4.9) 


-ia  ’y  -j-ia  'y 

0(-y  ,z)  -  Ie  °  °  cosh  a  z  +  Re  °  °  cosh  a  z  ,  0  s  2  s  h 

(4.10) 

Now  the  contribution  to  Eq .  (4.8)  from  z  =  0  vanishes  by  virtue  of  Eq.  (4.2). 
Also,  Eq.  (4.3)  implies  that  the  contribution  from  z=H  ,  -  ®  <  y  <  0 
vanishes.  Thus  Eq.  (4.8)  becomes 


,ra  J*  0  ^(y>.H?dy  +  1m  J(fi  dz-lm  §|)  dz  =  0  (4.11) 

O  O  '  ’  'ft  rt  ft  In 


The  first  term  presents  the  most  difficulty,  so  this  will  be  left  to  the 
last.  Substituting  the  asymptotic  form  for  0  into  the  second  Integral 
gives 


lm  j  -ib^coshPb^K  dz  =  2rj~(2boH  +  sinh  2bQH) 
o  Wo 

Similarly  for  the  final  term  in  (4. IT),  we  find 


h  / 
-im  J  -ia^: 


-ia'y  +ia'y\/_  +ia'y  _  -ia'y\ 


oKl  e  °‘°+  Re  0  °)(l  e  00  -  R  e  00jcoshaaoH  dz 


=  a  (2a  H  +  sinh  2aQH)  (|l|s  -  J  R  f S ) 
o 
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The  first  term  in  Eq.  (4.11)  may  be  written 


.  M  r°(5  ■  a,  30  3 {5 (y ,H)  . 

ln' k  J  (rr  - k )  sf  sr^dy 

O  Oy 


where  Eq.  (4.4)  has  been  used,  and 


.  30  30  . 

lmJ  eTaIdy  =  0 


It  is  possible  to  integrate  this  expression  by  parts  twice  to  obtain  the 
following  result 


j°(5f_.k>f  «4 1  dy 

\  3  oz‘  si:  9 

o  oy 


-  Im  K 


H[3y3W;H 


«[’ 


i  ni  .3,30  V  3  ,00, 

-  Im  —  |vk  (-T — ; _  -r —  (-X —  ;  , 
K  Oz  H  oy  'oz  ri 


(2~v) k2 


si  f|&  T 

S7''d2;H  ,  a  '5z^[. 

'  oy  j< 


where  all  the  integrated  terms  were  purely  real  and  hence  vanished  when 
the  imaginary  part  was  taken. 

Now  the  first  term  vanishes  at  y  =  0  from  condition  (4.6)  and  the 
second  term  in  brackets  is  seen  to  be  purely  real  from  Eq.  (4.5)  and  so 
it  also  vanishes.  The  value  of  this  expression  at  the  top  limit  yQ  may 
be  determined,  after  some  algebra,  to  be 


-f.b’b  3  sinh  2b  H 

o  o _ o 

l-L  +  Mb  4 
o 

Combining  the  contributions  from  each  of  the  terms  in  Eq .  (4.8),  we  find 
that 
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Mb  'b J3  sTnh  2b  H  b  ' 

- 2_e - S_  .  _2_  ^2b  H  +  s  inh  2b  H> 

,  ,  4  4b  x  o  o 

l-L  +  Mb  o 


+  {2aQH  +  s inh  2aQH)  (jl|E  -  |r|b)  =  0 
o 


This  may  be  written 


2b  H+sinh2b  H(l+ 
o  o 


4Mb  4 
o 

1-L+Mb  ‘ 


-)1  — -  -2-  (2a  H 

J  III2  % 


+  sinh  2a  H)(l  -  -f—M 

°  III2 

(4.12) 


Now  from  Eqs.  (3-43)  and  (3-44), 


*2  .1*1! 

Ill2 


and 


r*  =  la 


b  E  sinh3  b  H  , 
o  i 


a  sinh3  a  H  |lj‘ 
o  o  ' 


b  sinh  2b  K 
o  o _ 


a  sinh  2a  H  (1-L+Mb  *)  |l|; 

o  o  o  1  1 


so  that  in  terms  of  ^  and  ,  Eq.  (4.12)  may  be  written 


[>rB  +  /€2  =  1 


(4.13) 


where 


a  3  b‘  [(1-L+Mb  4)2b  H/sinh2b  H  +  (I -L+5Mb  4  )  1 
_  o  o  L _ o  o _ o  o  J 

3  2a  H / s inh2a  H  +  1 1 

o  °  J 


b  2 
o 


(4.14) 
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The  Case  M  =  0 


In  this  case 


where 


r  2b  H  -1 

si  nh  2b  H  1 
L _ o 

r  2a  H 

- 2 -  +  1 

si  n  h  Za  H 
L  O 


K  =  a  tanh  a  H 
o  o 


— —  =  b  tanh  b  H 
l-L  o  o 


so  that 


a  2b  1  [O-L)2  b  aH-kEH+K(l-L)l 

o  o  JL _ o _ J 

a  'b  2  1"  a  aH  -  K2H  +  k! 

o  o  L  O  J 


which  agrees  with  the  value  obtained  by  Weitz 


and  Keller 


3 


(05) 


(4.16) 


5.  SHALLOW-WATER  APPROXIMATION 
(a)  Formulation  and  Solution 

It  is  clear  that  the  full  solution  of  the  linearized  problem  presents 
overwhelming  numerical  difficulties  and  it  is  desirable  to  examine  the 
solution  on  the  basis  of  the  shallow-water  approximation.  Instead  of  con¬ 
sidering  the  limit  of  the  "finite  depth"  solution,  it  is  desirable,  for 
the  sake  of  clarification,  (and  certainly  easier)  to  rederive  the  solution 
from  the  original  approximate  equations  based  on  the  shallow-water  theory. 

We  return  to  the  Eqs.  (l.l),  (1.3),  and  (1.8)  satisfied  by  the 
velocity  potential  $(x,y,z,t)  . 

$+§+}=  0  ,  0<z<H  (5  1) 

xx  yy  zz  3  *  ^ 
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$  =0  ,  z  =  0 

z 


(5.2) 


4tt  +  g$z  =  0  ,  Z  =  H,  -  m<y<0,  -  co  <  X  <  ®  (5.3) 


#  +  =.5/^  -  —  h  $  ,  z-H  0<y< 

tt  3  z  p  z  p  ttz  ’  ’  ' 


_  a>  <  X  <  c  (5-4) 


We  have  seen  that  additional  conditions  at  the  leading  edge  of  the  ice-field 
are  given  by  the  equations 


_a_ 

dy 


(H*P>*(i-v) 

5y  5x§ 


z  =  H  ,  y  =  CM- 


£s  +  V 

ay2 


dxs 


=  0 


z  =  H  ,  y 


0+ 


(5.5) 


(5.6) 


where  §(x,y,t)  is  the  elevation  of  the  ice.  These  equations  ensure  that 
no  energy  is  put  in  or  taken  out  of  the  ice  at  the  leading  edge.  As  before, 
it  is  convenient  to  write 


5(x,y,z,t)  =  Re  -[0(y,z)e,kx  *  (5.7) 


so  that  0(y,z)  satisifies 


— ®  +  ~ — ^  -  k2  0  =  0  ,  0<z<H,-®<y<0  (5 .8) 

ay2  Sz2 


f?=0  ,  2  =  0  ,  -«<y<»  (5.9) 
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K0=  ||  t  Z=H  ,  -*  <  y  <  0  (5.10) 

K0  =  (1-L)  M  (^-  -  k2)  Sa  ,  z  =  H,0<y<»  (5-H) 

Sz  ay*  dz 

~  (If)  -  vks  (ff)  -  0  ,  z  =  H  ,  y  =  Of  (5.12) 

ay*  az  ez 

(|f)  -  (2-v)k2  f-  (|f)  =  0  ,  z  =  H  ,  y  =  0+  (5. 13) 

.3  az  ay  oz 

In  employing  the  shallow-water  approximation  the  potential  0  is  ex¬ 
panded  in  powers  of  z  ,  and  it  is  found  that 


0(y,z)  =  t(y)  -  jj;  z * (  ^yy  “  kS^  (5.14) 

satisfies  Eqs.  (5.8)  and  (5-9)  up  to  order  zJ  .  A  precise  derivation 

14 

of  the  theory  is  given  by  John  who  concludes  that  the  theory  is  valid 
provided  the  depth  of  the  water  is  small  compared  to  a  wavelength  and  to 
the  minimum  radius  of  curvature  of  any  immersed  body.  If  the  expression 
(5-14)  for  0(y,z)  is  substituted  into  (5.10)  and  (5-ll),  it  is  found 
that  'l'(y)  satisfies 

(i_  _  k2)<|i  +  =  0  ,  —  <  y  <  0  (5-15) 

Sy 

M(^—  -  k2)  t+  (l-L)  (■* —  -  k2)  t+  =  0  ,  Q<y<“  (5. 16) 

ay8  ays 


where 
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2  p;  hHK  c.  3 

K2  =  K/H  -  ,  L  =  -l--2-  ,  M  -  - ^ - 

°  9H  P  12(1-vE)d9 


in  terms  of  ^ (y)  the  conditions  (5-12)  and  (5.13)  become 


(-4  -  vk2)(-^-  -  k2)  t(y)  =  0  ,  y  =  04- 


(5.17) 


^  (-^  -  (2-v)k2)  -  k2)  <Ky)  =  0  ,  y  =  Of 

dy  5y2  Sy2 


(5.18) 


In  addition,  from  the  shallow-water  approximation. 


t(y)  and  ^  are  continuous 
dy  at  y=C 


(5.19) 


iK  'y  -iK  'y 

Now  the  Eq.  (5.15)  has  the  solution  t(y)  =  e  +  Re  ° 


where 


K  1  =  (K  -  k2 )  and  the  positive  square  root  is  taken,  and  K  >  k  for 

o  o  r  o 

a  progressive  wave.  The  form  of  the  solution  is  such  as  to  represent  an 
incident  wave  of  prescribed  amplitude,  and  a  reflected  wave  of  (complex) 

|  |^p  *  y 

amplitude  R  .  Eq.  (5.16)  has  a  solution  of  the  form  e  where 


k  ' 
n 


<kn  -  k  ^ 


(5-20) 


and  the  k  (n  =  0, I,  ...  5)  are  the  roots  of  the  equation 


Mk  6  +  ( 1-L) k  2  -  K  2 
n  no 


(5.2i) 


This  equation  is  a  cubic  equation  in  kn  having  solutions 


k  =  A  +  B 
o 


(5.22) 

(Cont 'd) 
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kj3  =  eA  +  ezB 
kg3  =  e®A  +  eB 


(5.22) 


where 


and 


2  1/3 
W 

V  27HK  4  / 

o  J 

a  1/3 

r  V8~ 

(V) 

V  27MK  -  ) 

o  J 

e  =  e2n,/3  =  (-1  +  -J3*  i )  /2 


(5.23) 


(5.24) 


In  all  eases  of  practical  interest 
A  +  B  >  0  . 


L  <  1 


and 


MJb 01 

27HKo4 


<  1 


so  that 


Thus  k^2  is  real  and  positive  while  Ic,2  =  kj2  . 


•  5  -> .  2 

The  square  root  in  the  expression  k  '  =  (k  -  k  ;  is  chosen  such 

n  n 

that  k  1  =  k  when  k  =  0  .  Now  for  y  >  0  a  solution  is  required 
n  n  _  _  ik0'y 

representing  a  progressive  wave.  Clearly  such  a  solution  is  e 
Other  solutions  are  permitted  provided  they  are  bounded  as  y  —  +  . 

Energy  considerations  exclude  the  possibility  of  unbounded  solutions.  The 
only  possibilities  are  those  roots  which  have  positive  imaginary  parts, 

so  that  e  n  ^  decays  exponentially  for  increasing  y  .  There  are  two 
such  roots.  ^Thus  let  k^  =  +  (kx2)2  have  a  positive  imaginary  part.  Then 
l<2  =  -  (kg2)‘  also  has  a  positive  imaginary  part,  since  Ic^  =  -  kj  . 

Hence  the  solution  for  y  >  C  may  be  written 


Ik  1  y  ik^'y  ik'y 

t(y)  =  A^e  °  +  AJ  e  +  A^e  2 


(5-25) 
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The  unknowns  are  the  complex  constants  A^,  Ax  ,  and  R  ,  and 
there  are  just  four  conditions  to  be  satisfied  by  t(y)  .  Since  i|r  and 
dt/dy  are  continuous  at  y  =  0  ,  then  the  following  equations  must  hold. 

Aq  +  \  +  A2  =  1  +  R  (5-26) 

ko'Ao  +  ki'\  +  =  Ko-(l  -  R)  (5.27) 

Also,  conditions  (5-17)  and  (5-18)  give  the  equations 


2 

^  Fk.3  -  (l-v)k2]  A.  =  0  (5.28) 

i=0 

2 

^k.'  k.s  Fk.s  +  ( 1  -v) ks ]  A.  =  0  (5-29) 

i=0 

which  together  with  Eqs .  (5.26)  and  (5.27)  are  sufficient  to  determine 
the  unknowns.  The  problem  is  solved  once  the  constants  Aq,  A1 ,  A2  and  R 
have  been  computed,  and  expressions  may  be  derived  for  physical  quantities 
of  interest. 


(b)  Reflection  and  Transmission  Coefficients 

The  transmission  coefficient  is  the  ratio  of  the  amplitude  of  the 
transmitted  wave  to  the  amplitude  of  the  incident  wave  at  infinity. 

Now  for  large  positive  y  , 

ik  'y 

(.(y)~Aoe  °  (5.30) 

since  the  other  two  terms  in  Ea.  (5.25)  decay  exponentially  with  increas¬ 
ing  y  .  The  elevation  S(x,y,t)  satisfies  the  Eq.  (1.6),  namely 
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S|  _  _  a« 

dt  dz 


2  =  H 


and  in  terms  of  i|i(y)  the  elevation  may  be  written 


Re 


'H  ,  2a 
—  k  A  e 
oj  o  o 


ikx 


+  i!<oy- 


and  the  correspond ing  amplitude  is 


(5-31) 


(5.32) 


The  incident  wave  is  given  ♦ 
cident  wave  is 


e 


so  that  the  elevation  of  the  in- 


?(*..  y>t) 


-  Re 


4"—  K' 

1(0  i 


ikx  + 


-  UaU 


(5-33) 


with  amplitude  .  Note  that  this  term  contains  the  dimensional  unit 

amplitude  of  the  incident  wave  and  hence  it  has  the  correct  dimensions. 
Thus  the  transmission  coefficient  if  is  given  by 


k  2 

<  =  —  I A  1  (5.34) 

K  3  0 


and  in  a  similar  manner,  the  reflection  coefficient 


H 


is  given  by 


£  =  |R| 


(5.35) 
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(c)  Pressure  on  the  Bottom,  z  -  0 

From  Eq.  (1-5)  the  pressure  on  the  bottom  in  excess  of  atmospheric 
pressure  pQ  is 


p  (g^-)  +  pgH 


z=0 


In  terms  of  the  shallow  water  potential  My)  >  we  have 


(5-36) 


pfx^yjOjt)  =  p  Re  •J-ittM(y)elkx  IU,tj-  +  pgH  (5.37) 

If  the  local  effects  represented  by  the  exponentially  decaying  terms  in¬ 
volving  Ax  and  in  Eq.  (5-25)  are  ignored,  then  the  amplitude  of 

the  pressure  fluctuation  on  the  bottom  under  the  ice  is 

<  ikx+ik'y-iuJt) 

p(x,y,0,t)  -  pgH  =  Re j  Pe  0  |  (5.38) 

where  jp[  =  piu  | | 

It  is  convenient  to  non-dlmens iona 1 ize  this  in  terms  of  the  amplitude  of 
the  incident  wave.  Thus  we  define 


r- 


9HK 


J  A 


s  i  nee  K  2  =  <“3 /gH  . 

o 

Hence  (P  is  a  measure  of  the  absolute  value  of  the  amplitude  of  the  pres¬ 
sure  fluctuation  on  the  bottom  z  =  0  ,  under  the  ice,  compared  to  the 
amplitude  of  the  pressure  fluctuation  on  the  bottom  under  the  free  surface. 
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(d)  Relation  Between  ^  and  T  for  Shallow  Water 

The  Eqs .  (5-26),  (5-27),  (5-28),  (5-29)  which  determine  the  coeffi¬ 
cients  R  and  A  conceal  the  relation  existing  between  and  l7~  . 
o 

This  is  most  easily  derived  by  taking  the  limit  for  small  H  of  the  rela- 
tion  which  has  been  derived  on  the  basis  of  the  full  linear  equations. 

From  Eg.  (^.13)  we  have 


o  -r- 


-fLz  =  l 


where 


D  =  —  —7 

b2  ao 
o 


(1-L+Mb  4  )  2b  H/s i nh2b  H  +  (1-L+5H  b  4 )  j 
o  O  O  _ O  J 


^2aQH/s ] nh2a^H  +  lj 


so 


2b  H 

H  -  0  >  TT^Th  -  1  +  0(boH)S 


a  4  b‘ 

D  -  —  —  (1-L  +  3Mb  4) 

b  2  a‘  0 

o  o 

where  b  is  the  real  positive  roots  of  the  equation 
o 


K  =  b  (l-L  +  Mb  4 )H 

o  o' 


since  tanh  b  H  ~  b  H(1  +  0(b  h)  * 

00  o 


This  equation  may  be  written 


Mbo6  +  0-L)do3  =  -here  ^=^=Ko 
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Clearly  b  —  k  in  the  shallow-water  notation  and  b  1  -*  k  1 .  On  the 
o  o  o  o 

same  basis 

K  =  a  SH  so  that  a  -  K  and  a  '  —  K  ' 
o  o  o  o  o 


Thus 


C 1 — L  +  3Mk  4) 
o 


and  for  shallow  water 


K  2k  ' 


k  EK  ' 
o  o 


( 1  -L  +  3«k  4)  “f2  +  fl* 


1 


If  we  introduce  the  incident  angle  8  and  a  transmitted  angle  by 

writ ing 


k  =  K  s  in  9  -  k  sin  B. 
o  o  T 

then  the  above  formula  becomes 


K  cos  9. 
o 


k  cos 
o 


g1  (1-L  +  3Mko4)rS  +  iC'  =  1 


(5-39) 


This  relation  will  provide  a  check  on  the  computed  values  of  £  and  . 


(e)  The  Critical  Angle  in  Shallow  Water 


As  in  the  case  of  finite  depth  of  water,  the  constant  k  in  the 

j 

x-variation  e  determines  the  angle  at  which  the  incident  wave  enters 

tiie  ice-field.  Thus  if  we  write  k  =  K  sin  8  ,  then  the  incident  wave  is 

o 

of  the  form 


e 


iK  y  cos  6  +  iK  x  s in  6 
o  o 
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so  that  the  incident  wave  makes  an  angle  3  with  the  normal  to  the  leading 
edge  of  the  ice-field,  as  shown  in  Fig.  2. 

The  reflected  wave  is  of  the  form 


-iK  y  cos  6  +  IK  x  sin  6 
o  o 


indicating  that  the  angle  of  incidence  is  equal  to  the  angle  of  reflection. 
The  transmitted  wave  has  the  form 


i(ks-K2sinsS)2y+iKxsinS  ik  cos  By  +  ik  sin  @Tx 

oo  o  _  o  T  o  t 


where  k  satisfies 
o 


n  is. 

i<(  6  .  i  \  i  S  i,  3  f  sin  o  o 

Mk  +  (l-L)k  =  K  ,  and  -  7 T" 

o  o  o  sin  9-  K 

T  o 


(5-40) 


Clearly  when  ®  =  0  ,  sin  0  =  0  ,  and  since  kQ  is  real,  there  will 

always  exist  an  undamped  transmitted  wave.  (Assuming  L  <  1  .)  But  if 

for  particular  values  of  M  and  L  ,  K  >  k  then  there  will  exist  a 

o  o 

critical  value  of  6  at  which  k  3  =  K  3  sin3©  and  above  which 
x  o  o 

(k^3  -  KqS  sin3  ©)2  is  pure  imaginary,  thus  producing  an  exponentially 

damped  transmitted  wave.  This  value  of  9  is  given  by  9  .  =  sin  1  (k  /K  ). 

v  a  7  crrt  o  o 

An  incident  wave  approaching  the  ice  at  an  incident  angle  greater  than 
0  will  be  totally  reflected  by  the  ice  as  shown  in  Fig.  3. 

On  the  other  hand,  if  k  >  K  then  (k  3  -  K  a  sin3  B)  is  always 

oo  oo 

positive,  so  that  waves  approaching  at  any  angle  will  penetrate  the  ice. 

In  this  case  the  transmitted  wave  is  bent  towards  the  normal  so  that  the 

ice  may  be  regarded  as  being  "denser"  than  the  water.  When  K  «=  k  , 

o  o 

which,  from  Eq.  (5.4c),  occurs  when  MKq4  =  L  ,  the  transmitted  wave  pro¬ 
ceeds  at  the  same  angle  as  the  incident  wave.  The  wavelength  at  which 
this  occurs  is  denoted  by  ^-crjt  and  hence  any  wave  satisfying  A  >  A.^. 
penetrates  the  ice  regardless  of  the  incidence  angle. 
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6.  DISCUSSION  OF  RESULTS 


(a)  Description  of  Procedure 


In  obtaining  the  numerical  results,  the  following  values,  taken  from 


Robin,  were  used 


Young’s  modulus  for  ice  E  =  SxlO1  dyn/cm 

Poisson’s  ratio  for  ice  v  =  0.3 


Density  of  ice 


p.  -  0.92  gm/cm3 


Density  of  sea -water  p  =  1.025  gm/cm 


„  2tt  ,  ^i  /2n%2  35-4Hh 

Ko  =  T  '  50  that  L  =  T  (T>  = 


H= - - -  =  4.54  x  106h3 

12(l-v3)pg 


Robin'  has  observed  ice  thicknesses  of  about  1.5  metres  so  that  values  of 
h  =  0.75,  1.5,  3-0  metres  were  considered  together  with  water  depths  of 
H  =  10,  20  metres.  The  incidence  angle  9  was  allowed  to  vary  in  steps  of 
15°  from  0  to  60°,  and  various  wavelengths  of  the  incident  wave,  up  to 
700  metres,  were  considered. 

In  each  of  these  cases,  the  Eqs.  (5.26),  (5.27),  (5.28),  and  (5.29) 
with  coefficients  determined  from  Eqs.  (5.22),  (5.23),  and  (5.24)  were 
solved  and  values  were  obtained  for  '7',  ,  and  fi  . 

In  addition,  the  critical  angle  ®crjt  above  which  an  incident  wave 
is  completely  reflected  and  the  critical  wavelength,  A.  .  above  which 
an  incident  wave  at  any  angle  penetrates  the  ice,  were  computed  in  each 


As  a  check  on  the  numerical  work,  the  expression 


K  cos  9_(  l-L+3Mk  4)f3 
o  I _ o _ 

k  cos  0 
o 
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was  evaluated.  This  should,  of  course,  from  Eq.  (5.39)  be  identically 

zero.  The  fact  that  this  expression  was,  in  all  cases,  negligibly  small 

provided  a  check  on  the  numerical  accuracy  as  well  as  a  verification  of 

the  shallow-water  approximation  since  the  relation  between  °T  and  ^ 

was  derived  by  taking  the  limit  for  small  H  of  the  corresponding  result 

based  on  the  exact  1  inear  theory.  Presumably  this  same  relation  could  be 

derived  directly  from  the  analytic  solution  of  the  equations  determining 

A  and  R  . 
o 

The  results  are  shown  in  Figs.  5  through  15- 
(b)  The  Critical  Angle 

In  Fig.  5,  the  critical  incident  angle  is  shown  as  a  function  of 
incident  wavelength.  The  critical  angle  determines  whether  or  not  a  given 
incident  wave  at  a  given  incident  angle  is  able  to  penetrate  the  ice  (see 
Figs.  2  and  3).  Thus  a  wave  approaching  the  ice-field  at  an  angle  9 
greater  than  the  critical  angle  (as  determined  from  Fig.  5)  is  completely 
reflected  by  the  ice  field.  If  the  incident  anqle  is  less  than  0 
then  penetration  of  the  ice  field  occurs  and  an  undamped  transmitted  wave 
travels  through  the  ice  field.  Incident  wavelengths  greater  than  4  . 

corresponding  to  8  .  =  will  always  penetrate  the  ice,  regardless  of 

the  incident  angle  (see  Fig.  4). 

(c)  Reflection  and  Transmission  Coefficients 

Figures  6  to  10  show  the  variation  of  and  *T  with  the  incident 
wavelength  X  for  ice  thicknesses  of  0.75,  1-5,  and  3.0  metres,  incident 
angles  ranging  from  0°  to  60°  in  steps  of  15°,  and  water  depths  of  10  and 
20  metres. 

In  Fig.  6  the  curves  for  7”  have  been  extrapolated  so  as  to  pass 

through  the  origin.  This  is  not  strictly  accurate  since  at  X  =  0  ,  L 

is  infinite,  and  it  has  been  assumed  that  L  <  1  .  There  does  in  fact 

exist  a  cut-off  wavelength  (unrelated  to  the  incident  angle)  which  occurs 

whenever  Eq .  (5.21)  fails  to  have  a  real  root.  This  occurs  when  ^jj-  >  1 

L  /rm  * 

and  then  all  possible  solutions  in  the  ice-field  decay  exponentially  0 
with  distance.  This  restriction  is  of  little  physical  significance  since. 
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for  example,  with  h  =  1.5  metres,  H  =  10  metres,  no  propagation  is  possible 
for  A.  <  10~*  cms,  approximately.  In  other  words,  this  mathematical  phenom¬ 
enon  only  occurs  at  wavelengths  too  small  to  be  physically  significant. 

Figures  6  to  10  illustrate  the  importance  of  ice  thickness  in  deter¬ 
mining  reflection  and  transmission.  Thus,  for  instance,  from  Fig.  6, 
doubling  the  ice  thickness  from  1.5  to  3  metres  causes  a  drop  of  over  25% 
in  the  wave  amplitude  of  the  transmitted  wave  corresponding  to  an  incident 
wavelength  of  250  metres.  Since  the  wave  energy  is  proportional  to  the 
square  of  the  wave  amplitude,  this  means  almost  a  50%  reduction  in  energy 
caused  by  doubling  the  Ice  thickness. 

In  contrast,  the  water  depth  has  little  effect  on  the  reflection  and 
transmission  coefficients  except  for  the  smaller  wavelengths. 

As  would  be  expected,  the  ability  of  the  ice  to  reflect  the  incident 
wave  is  reduced  as  the  incident  wavelength  increases,  and  for  A  >  500  metres, 
the  transmitted  wave  height  is  at  least  95%  of  the  incident  wave  height. 

In  all  cases,  the  amount  of  reflection  was  very  small  except  for  the  lower 
wavelengths.  This  agrees  with  the  conclusions  of  Stoker,^  who  considered 
the  two-dimensional  problem  (6  =  0°). 

Figures  6  to  10  furthermore  indicate  the  effect  of  the  angle  of 
incidence  upon  the  reflection  and  transmission  coefficients.  For  a  fixed 
wavelength  there  appears  to  be  a  gradual  increase  in  transmission  with 
increasing  incident  angle.  For  a  given  incident  angle  '7~  decreases  until 
the  cut-off  wavelength  is  reached,  at  which  "T  drops  to  zero  and  complete 
reflection  occurs.  For  larger  incident  angle  a  distortion  of  the  trans¬ 
mission  curves  takes  place  so  that,  for  example,  when  9  =  60°,  the  trans¬ 
mission  curve  for  ice  of  3  metres  thickness  has  a  minimum  at  a  wavelength 
of  about  320  metres. 

(d)  Pressure  Amplitude  on  the  Bottom  Under  the  Ice 

Figures  11  to  15  show  the  variation  in  non-dimensional  pressure 
fluctuations  on  the  bottom  under  the  ice  normalized  with  respect  to  the 
pressure  on  the  bottom  in  the  absence  of  ice,  in  terms  of  the  incident 
wavelength  \  ,  ice  thickness  h  ,  water  depth  H  ,  and  incident  ancle  9  . 
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In  all  cases  the  pressure  amplitude  ft  tends  to  unity  as  4  in¬ 
creases,  which  is  to  be  expected.  For  6  =  0°,  the  pressure  is  generally 
lower  than  the  corresponding  pressure  amplitude  on  the  bottom  in  the  ab¬ 
sence  of  ice,  the  pressure  drop  being  greater  for  greater  ice  thickness. 

For  6  =  15°,  considerable  changes  in  the  pressure  curves  take  place  and 
for  the  lower  incident  wavelengths  P  is  actually  greater  than  unity  for 
an  ice  thickness  of  3  metres.  This  trend  increases  for  larger  values  of 
8  so  that  for  9  =  60°  the  pressure  jumps  to  over  1.5  the  free  surface 
bottom  pressure  in  the  most  extreme  cases.  As  in  the  case  of  the  trans¬ 
mission  coefficients,  for  each  non-zero  angle  there  is  a  cut-off  wavelength 
corresponding  to  0  .  at  which  the  pressure  drops  abruptly  to  zero, 

indicating  complete  reflection  of  the  incident  wave. 

(e)  Comparison  of  Results  With  Observation 

Robin^  has  concluded  that  no  effective  transmission  of  wave  energy 
occurs  for  wavelengths  less  than  200  metres  in  fields  of  large  floes, 
whereas  major  energy  changes  occurred  in  those  waves  having  a  period  of 
16  seconds,  or  wavelengths  of  400  metres. 

It  is  difficult  to  compare  observations  in  deep  water  with  a  theory 
based  on  the  shallow  water  approximation.  Thus  Figs.  6  to  10  indicate 
that  energy  transmission  through  the  ice  occurs  at  much  lower  wavelengths 
than  ou^erved  by  Robing  For  instance,  when  A.  =  100  metres,  h  =  3.0m, 
it  =  0.3  indicating  that  over  90%  of  the  incident  energy  has  been  trans¬ 
mitted  into  the  ice.  This  discrepancy  is  almost  certainly  due  to  using 
the  shallow-water  approximation.  If  the  numerical  difficulties  involved 
in  the  solution  based  on  the  full  linearized  theory  could  be  surmounted, 
it  is  felt  that  a  much  closer  comparison  of  theory  and  observation  would 
result. 

The  experimental  counterpart  of  this  wave-ice  investigation  is  re¬ 
ported  in  Ref.  15. 
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B.  MODEL  ID.  AN  ICE-FIELD  HAVING  SURFACE  TENSION 
1.  FORMULATION  AND  SOLUTION 

In  this  model  it  is  assumed  that  each  element  of  the  ice  sheet  is 
subjected  to  a  force  arising  from  a  surface  tension  force  in  the  sheet. 
Then  it  may  be  shown  that 


P  -  Pc  =  -  Tsv£§  +  Pjh§tt  (6.1) 

where  Ts  =  surface  tension  force. 

The  theory  proceeds  exactly  as  for  Model  II  except  that  the  free 
surface  condition  satisfied  by  iJj  becomes 


°y 


z  =  H  ,  0  <  y  < 


where  S  =  T  /pg  . 


The  characterist ic  equation  is 


(6.2) 


K  cos  b^H  +  bn  jl-L-Sbn2j  sin  bnH  =  0  (6.3) 

which  has  an  infinite  sequence  of  real  roots  ±bn(n  =  1,2,  ...)  and  two 

pure  imaginary  roots  ±ibQ  .for  K.S  >  0  ,  L  <  1  .  Thus  propagation  always 

occurs  in  the  ice,  if  | f  >  k  .  Also,  complete  reflection  of  a  plane 

incident  wave  will  always  occur  for  some  angle  of  incidence  whenever 

L  <Sa  2  . 
o 

The  transform  of  condition  (6.2)  is  obtained  by  integrating  by  parts; 

thus 

KY+  =  Y+'(l-L  +  SV2)  +  SCa^iOaJ  (a  e  D+)  (6.4) 
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with  the  same  notation  as  before.  Similarly,  for  y  <  0 

K*_  =  Y_'  (or  e  D_)  (6.5) 

Following  the  same  procedure  as  in  Model  1 1 ,  we  arrive  at  the  equation 


K+(a)  j(SYa-L)*;  -S(a3-iffl)J  -1A 


.  (K+(Qf)-K+(-by 


QfH-b1 

O 


~M«)  j  (SY*-L)'l!:-S(as-iffal)J  -iA 


K  <«>-K+<-by 


«+b' 

o 


=  j(») 


(««D)  (6.6) 


where 

KcoshvM-v(l~L+S,Ya)sInhYH  K+(a) 

K  cos  hyH-Ys  i  n  hry H  =  K_  (a) 

and  K^C®)  is  regular  and  non-zero  in  D±  ,  respectively. 

The  left-hand  side  of  (6.6)  is  regular  everywhere  in  D+  while  the  right- 
hand  side  is  regular  everywhere  :n  D  The  two  sides  are  equal  in  the 
strip  D  .  This  defines  a  function  J(a)  regular  in  the  whole  a-plane. 

It  may  be  shown  that 

K+(o;)  =  0 (or)  lorl  -  «  in  D+ 

K_(or)  -  0(a-*)  lorl  — «  in  D 

Furthermore  T  1  =  o(  l),  1 0f|  —  “  in  ,  as  in  Model  II. 

But  the  left-hand  side  of  (6.6)  is  just 

=  •>(<*)  (6.7) 


K+(«)  {k*+  -  ¥+'|  - 


K+(«)  - 


K+(- 


b  ') 


+  b  ' 
o 
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and  so  J(a)  =  o(a)  as  |o'|  “  in  D+  . 

Similarly  from  the  right-hand  side  of  (6.6) 

J(a)  =  o(»)  ,  as  JcrJ  -  “  in  0 

Thus  from  an  extension  of  Liouvi lie's  theorem 

J(of)  =  C  ,  a  constant 

We  find  from  Eg.  (6.6),  that 


^(ar.H)  = 


C  - 


iAK+(- 


b  ') 


o  +  b  ' 
o 


_ 1 _ 

(K  -  ytanhyH)K+(cs') 


and  so 


*«■«>•  wL  *<»■">£§£  »-ioy« 


where  the  path  of  integration  is  in  D  . 

We  see  that  the  solution  depends  upon  a  constant  C  .  An  additional 
condition  is  required  to  determine  C  . 

The  required  condition  is  obtained  by  applying  Green's  theorem  to 
the  function  0(y,z)  and  its  complex  conjugate  0  in  the  rectangle  formed 


by  the  lines  z  =  0,z=H,y  =  ±yc 
that  0(y,z)  satisfies 


Note,  by  analogy  with  Model  II, 


(as/Sy2  +  62/bz2)0(y,z)  -  0 


(6-8) 


KJ5  = 


50 

37 


,  z=H,  ~®  <  y  <  0 


(6.9) 
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where  K  =  — 
9 


S 


T 


s 

PS 


Tb'y  . 

0(y,z)  =  T  e  °  coshb  z  +  0(e  ^),  0  ^  z^H,  y>0  (6. 1 1) 


Ji(y,z)  =  I  e'a°^coshaoz  +  R  e  °ycoshaQz  +  0(e^y) 

0^z^H,y<0 


(6.12) 


where 


i  4- 

a1  *=  (a  2  -k*j  ,  b*  =  (b  8  -k2) 
o  '  o  *  o  '  o  1 


and 

K  cos ha  H  =  a  sinha  H 
o  c  o 

K  coshb  H  =  b  (l-L+Sb  2)  sinhb  H 
o  o  o  o 

The  method  is  identical  to  that  used  in  Model  II  to  determine  the 
relation  between  and  ‘T  .  In  addition,  it  provides  us  with  the 
required  conditions  for  the  potential  at  infinity  to  be  unique. 

Green's  theorem  reduces  to 


ImJ  tf^ds  (6.13) 

C 

where  the  line  integral  is  taken  over  the  boundary  of  the  rectangle.  The 
contributions  from  z  =  0  ,  and  from  z=H  ,  -“<y<0  vanish,  and  the 
contribution  from  z  =  H,  Q  <  y  <  °°  is 
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y  y 

■  r°  *  .  .  -s  r°  .2wa«i.  ,30. 

.  k  )(57)z=h(5_)z=h  dy 

-j  y=y- 

t=H  _n 

Jv=0 


=  T  ""[i 


o 

3_,30 

y 


'  o 
Jy=0+ 


-$b2sinhab  Hb'  ~ 

°  °_°+£  im<|_|_|0) 

K  dyoz  °z  z=H 
y=0+| 


K 


T|S 


The  contribut ions  to  Eq.  (6.13)  from  the  boundaries  y  =  ±  y  are 

o 

ao  (2a0H+sinh2c0H)(jl|s-|R|3  ) 

4a 

o 


and 


-b^  (2boH+sinh2b  H)  )t|3 
—  ,  respectively 

O 


Thus  combining  all  the  terms,  and  using  the  relation 


Kslnh2b  H 

sinh2b  K  =  — - 2 - - 

°  2bo0-L-,Sbo-) 


we  obtain 


,  b1  f 

hi*  s2-  [ 


2Sb  45 

2b  H+s i nh2b  H(]+  - — 


1 -L+Sb 


-) 


v"l  -  is  ,  ,3S0 

>J  +~  ,m(^^}z=H 

y=0+ 


a'  (2a  H+sinh2a  H)(|ljS -|R|a) 


(6.14) 


Now  from  Eq.  (6. l4) ,  if  |l|-  -  0  ,  so  that  no  energy  is  being 
supplied  to  the  system,  R  =  T  =  0  if  and  only  if 
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"  0  (6‘,5) 

y=0+ 

This  condition  must  therefore  be  satisfied  to  ensure  that  no  energy  is 
being  supplied  to  the  system  other  than  by  the  incident  wave.  Then  the 
following  relation  may  be  derived  between  the  reflection  coefficient 


t€2 


(6.16) 


and  the  transmission  coefficient 


3 

T 


b  2sinh2b  H 
o _ o  _ 

a  2sinh2a  H 
o  c 


III3 

III8 


(6.17) 


dT2  +  j?3  =  i 


(6.18) 


where 


ao  bo  [(1-L+Sfa02  )2boH/sinh2boH  +  (l-L+BSb^)] 

b  2  a'  f  2a  H/s inh2a  H  +  ll 

o  o  L  o  o  J 


(6.19) 


When  H  is  small,  this  reduces  to 


_  o  o  (i-L+2Sk0® ) 


(6.20) 


in  the  notation  of  the  shallow-water  approximation  used  in  Model  II. 
The  condition  (6.15)  may  be  satisfied  in  more  than  one  way  and  the 
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physical  unreality  of  the  model  does  not  provide  an  obvious  choice.  How¬ 
ever,  the  tension  existing  in  the  ice  sheet,  z  =  H,  y  >  0  suggests  there 
must  be  a  restraining  force  at  the  edge  of  the  sheet  y  =  0+.  For  instance, 
if  the  ice  is  fixed  at  that  point  in  some  way  such  that  the  elevation  is 

zero,  then,  (4®)  =0  and  condition  (6.15)  is  satisfied.  This  condition 

°z  z=H 
y=0+ 

is  sufficient  to  determine  the  constant  C  so  that  the  solution  is  uniquely 
determined. 

This  work  provides  an  extension  of  the  problem  considered  by  Keller 
and  Goldstein3  who  utilized  the  shallow-water  approximation.  The  authors 
apply  a  different  third  condition  at  the  edge  of  the  ice-fieid  in  order 
to  fully  determine  the  solution.  ihis  condition  which  is  essentially  the 
continuity  of  vertical  velocity  at  the  leading  edge,  is  not  sufficient  to 
ensure  that  no  energy  is  fed  into  the  system  other  than  from  the  incident 
wave,  and  also,  the  authors  do  not  discuss  the  instability  of  the  stretched 
membrane  at  the  edge  z  =  H,  y  =  Q+. 
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APPENDIX  A 


It  is  required  to  find  the  position  of  the  roots  of  the  equation 


Kcosb  H  +  b  -fl-L+Mb4}  sinb  H  =  0,  L  <  I,  M,  K  >  0 
n  n  l  nj  n 


(A-l) 


By  sketching  the  curves  tanb  H  and - p - —  ,  it  may  be  shown  that 

n  b  ]  i  -  L+Mb4  f 

n  k  nJ 

there  is  an  infinite  sequence  of  reai  zeros  -bn(n  =  such  that 

!  b  <  b  . !  ,  and  furthermore 
n‘  1  n+l 


bn  ~f  +  °<4^  aS  "  - 


(A-2) 


The  intersection  of  the  curves  tanh  b^H  and  ^  (1-L-tHb4) 
are  also  two  pure  imaginary  roots  -ib  - 


indicate  there 


The  number  and  position  of  the  complex  roots  of  (A-l)  are  deter¬ 
mined  by  considering  the  equivalent  equation 


A  +  z(l+Bz4)  tanz  =  0 


(A-3) 


for  complex  z(=  x  +  iy). 


Now 


z(l+Bz4) 


0  as  I  z 


whereas 
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i  i  /  S  i  n*  x  +  s  i  nh*  y\ 

1  tanz  |  -  [ - - - — 

\COS^X  +  S  I nh  y / 


-0 


provided  x  —  rrrr  and  y  —  0,  simultaneously- 

Equation  (A~3)  may  be  split  into  real  and  imaginary  parts  and  one 
form  of  the  resulting  equations  is 


~  Yt)  (dgX8  -  y3)  s  i  nh2y  _  _  sin2x  (A_4) 

(eix2  -  ^)(e2xa  -  /)  2y  2x 

A^sjj2h2jr  +  S"',^X)  =  2B{x^  -  y4)  (cosh2y  -  cos2x)  (A-5) 


where 


dt  ,s  =  1  -  V20  -  5/Bx4 

et  >8  =  5  -  yJlo^T/Bx* 


Now  as  x  —  nrr,  n  large,  y  -»  0,  the  left-hand  side  of  (A-4)  tends  to  I, 

while  the  right-hand  side  tends  to  zero-  Hence,  there  exists  no  solution 

to  (A-3)  for  x  >  x  for  some  x  =  x  (A,B)-  Also,  it  is  clear  from  (A-5) 
c  o  o 

that  x>  y  for  a  solution  to  exist-  Now  from  (A-3),  if  zo  is  a  solution 

with  o  <  erg  z  <  t  >  then  -z  ,  z  and  -z  are  also  solutions-  There  are 
4  o  o  o 

thus  a  finite  number  of  solutions  of  (A-2)  occurring  in  pairs  within  the 
region  R  bounded  by  the  lines  y  =  -x,  x  =  rmr  ,  for  n  large  enough- 
(See  sketch  on  the  following  page). 
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The  number  of  roots  lying  in  R  may  be  determined  by  the  method  of  the 
argument-  Thus,  if 


f(z)  =  A  +  2(1+Bz4)tanz 


— —  A  arg  f(z)  =  (no-  of  2eros)  -  (no-  of  poles) 
2n 


where  A  arg  f(z)  means  the  change  in  the  argument  of  f(z)  as  z  traverses 
the  boundary  of  R.  It  may  be  shown  that  ^  A  arg  f(z)  =  2  and  since  there 
are  n  real  zeros  and  n  poles  T n  R,  then  there  must  also  be  two  complex 
zeros.  Thus,  in  tbe  entire  plane  there  are  four  complex  zeros  of  f(z), 
such  that  x  >  y- 

Thus,  Eq.  (A-l)  has  complex  roots  ±c  ,  ±c  where  c  =  X  +  ip 

□  o  o 

and  \  >  p  >0- 


In  a  similar  manner,  it  may  be  shown  that  the  equation 


K  cosa  H+asinaH=0  (A-6) 

n  n  n  v 

has  an  infinite  sequence  of  real  zeros  ap  (n  =  i,2,---)  such  that 

a  <  la  ,  ,|,  and  furthermore 
n  n+ 1 
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a 

n 


~  T+0^) 


as 


n  -*  w 


(A-7) 


There  are  also  two  pure  imaginary  roots  -ia  -  It  may  be  shown  that  there 
are  no  complex  roots  by  considering  the  equation 

A  +  ztanz  =  0  for  complex  z  =  x+iy  (A-8) 


When  (A-8)  is  split  into  real  and  imaginary  parts,  it  is  found  that  one  of 
the  resulting  equations  may  be  written 


sinh2y  sin2x 
2y  +  2x 


(A— 9) 


Now 


s  inh2y 
2y 


whe re as 


sin2x 

2x 


<  1 


SO  that  (A-9)  has  no  solution  for  x,  y  ^  0. 
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APPENDIX  B 

It  is  required  to  determine  functions  K±(o0  such  that 

f,  H 

fM  “  K+H/K_(^ 

o' 

and  K±(it)  is  regular  and  non-zero  in  D_j_  . 

Consider  the  function 

fi  (<v)  =  K coshy H  -  y(l  -L+MyA)  s  i  nhyH 

Th is  has  zeros  y  =  ±b  ,  y  =  ±ib  (n  =  i ,2 , • - and  y  =  ±ic  ,  ±ic  , 
on  o  _  o 

Re  c  >  0-  Thus,  since  y2  =  ^  +  k?  .  n  =  ib 1 ,  ±ib',  ±ic: ,  ±ic‘  ,  where 
0  o  n  o  o 

Re  c^  >  k  since  0  <  arg  cq  <  7;  ■ 

Since  f,  (o')  is  an  entire  function  (a  function  whose  only  singular- 

1  I 

ities  are  poles),  it  may  be  written  as  an  infinite  product.  (Noble  ,  p .  1 5 ) 
Thus 


ft  W  =  K(1  -  ^4)(1  +  7^)0  +  -^r)  |  |  (1  +  t4)  (n  =  1,2,-.-) 

000  h=l  n 

where  y3  =  op  +  k2  . 

In  a  similar  manner,  the  function  f^ry)  =  KcoshyH  -  ksinhyH  may 
be  r i  t te -n  in  the  form 


f 

o 


(ry) 


=  Kd  -  —g 

O 

o 


)  ft  <1  *-£r) 


8! 


R-1313 


Once  the  singularities  of  fo(w),  (o')  are  revealed  in  this  manner,  it 
is  possible  to  define  K±(«)  in  the  following  way- 

Let 


K+(«>  «  0+fa)  -f 

o  o 


k2  x 2  in' 

{,  +  T^>  ‘  T 

O  O 


f,- 

n=I 


(i 


_  fo^+k3) 


K  (a)  = 


-  ■Iar£1f  jo*  fs-A-Hje*#.)1*  rj  1 1  j»*g>**f 

o  M  o  cm  o  o )  n-l  (  n  n 


1 

(i  +  ^V)  + 

a  a 
n  n 


Then 


^  _  KcoshyH  -  Y(l~L+MY4)sinhYH 
KcoshyH  -  ysinhyH 


and  K^(t)  is  regular  and  non- zero  in  D±>  respectively. 
It  is  noteworthy  that 


K  (o')K  (-«)  = 


.  (c^-f  ke 


b* 


and  that 


Mai> 


=  1 


The  behavior  of  K±(t*)  as  Ini  -  <=  in  may  be  determined  by  considering 
the  behavior  of  the  infinite  products  ss  Iwl  —  =  -  Now  it  may  be  shown  that 


rrrr  n  ,  K 

an=T  +  0(n)  aS  n 
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whereas 


Thus 


bn=T+0^> 


n=l 


I  it2  2 

f  n 


d+ J-)2  - 

a  a 


TT  <r> 

n=l 


CM-i(b®  +  ks) 


n=l 


®+i  (a®  +  k2  ) : 


Now 


Also 


'  +  0(n2)  so  that  "j  I"  (-jp)  =  constant  < 
n=  I  " 


n 

n=l 


"+?(b3  +  k2!3 
_ n 


»+  i  (a2  +  k3 ) : 


TT 

n-| 


{'  +  fnW} 


whe  re 


fnH  = 


<bn  +  k?>  -  (<  *  kS) 


(a2  +  k  )  -  i(v 


Clearly  f  M  ~  0  as  (  or  |  -  =  ,  rreD 


and  |f  (or)  1  <  &Bi£ 


so  that 
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,nm  TT  (i+f  («))  = 

|ry  |->  =>  n=  | 


FT  ';m  (i+f  w)  =  i 

n=l  hb  ® 


Hence  the  infinite  products  in  the  expression  for  K^o1)  are  0(1)  as 
bb  “»  »gD+-  A  similar  result  holds  for  the  infinite  products  appearing 
i  n  K_  (o')  ■  Thus 


K.+  (»)  =  0 (w3)  as  lo< I  -•  «  in  D+ 
K_(ry)  =  0 (o'  )  as  Inrl  —  es  in  D_ 
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TABLE  1.  FREQUENCY  AND  WAVELENGTH  BANDS  CORRESPONDING 
TO  INCIDENT  WAVES  WHICH  ARE  COMPLETELY  REFLECTED 
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FIG.  I.  STABILITY  DIAGRAM  FOR  SOLUTIONS  OF  MATHI EU'S 

DIFFERENTIAL  EQUATION  (REPRODUCED  FROM  McLACHLAN, 
THEORY  AND  APPLICATION  OF  MATHI  E  U  FUNCTIONS,  P.  40) 
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FIG.  2.  REFLECTED  AND  TRANSMITTED  WAVE  FOR  0<©CR|T 
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FIG.  3.  COMPLETE  REFLECTION,  NO  TRANSMITTED  WAVE  FOR  e>6CR|T 


FIG.  4.  ALWAYS  TRANSMISSION  FOR  X>X„D1T 
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FIG.  5.  RELATION  BETWEEN  CRITICAL  ANGLE  0CR|T  AND  INCIDENT  WAVE  LENGTH 
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